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Young's  modulus  of  the  test  Material 

Young's  modulus  in  compression  of  the  test  material 

Young's  modulus  in  tension  of  the  test  material 

The  probability  of  failure  of  a  component 

Shear  modulus  of  the  beam  material 

Moment  of  inertia  for  a  rectangular  beam  (I  «  bd3/12) 

Moment  of  inertia  for  a  rectangular  beam  with  45°  chamfered  corners 
(See  Appendix  F) 

Moment  of  inertia  for  a  rectangular  beam  with  round  corners  (See 
Appendix  F) 

Outer  span  length  for  a  four-point  and  a  three-point  loaded  beam 
Total  length  of  beam 

Weibull  slope  parameter,  the  "Weibull  Modulus"  associated  with  either 
volume  or  surface  sensitive  material 

General  moment  applied  to  beam 

Bending  moment  as  a  function  of  x  (See  Appendix  D  or  E) 

General  applied  force 

Forces  applied  to  a  beam  (See  Figure  1) 

Effective  surface  of  a  beam  in  bending 
Torque  associated  with  beam  twisting 

Estimated  torque  when  bottoming  of  the  load  fixture  occurs 
(See  Appendix  C) 

Volume  of  a  beam  in  bending 

Effective  volume  of  a  beam  in  bending  (See  Weibull  Analysis) 

Volume  of  a  three-point  loaded  beam  (VL  *  Lbd)  in  the  risk  of  rupture 
equation 

Half  the  distance  between  the  inner  span  and  outer  span  for  a  four- 
point  loaded  beam,  i.e.,  (L-i)/2  or  a»L/2  for  a  three-point  loaded 
beam  (note  a^  ■  a£  *  a) 

A  beam  dimension  (See  Figures  1  and  2) 

Beam  width  (See  Figure  1  or  5) 

Chamfer  of  a  comer  of  a  beam  vith  45°  chamfers  (See  Figure  5) 

Beam  depth  (See  Figure  1  or  5) 

Load  eccentricity  equal  to  (a* -a) 

Load  eccentricity  ratio  equal  to  (aj-a)/L 
Shift  of  neutral  axis  in  an  initially  curved  beam 
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Error  in  the  specimen  diemnsions  b  or  d 
Error  in  the  span  lengths  i  or  L 

Horizontal  shift  of  contact  and  load  points  due  to  beam  bending 
(See  Figure  4) 

A  numerical  factor  dependent  upon  b/d  (See  Appendix  C) 

Inner  span  length  of  a  four-point  loaded  beam  (See  Figure  1) 

Either  equal  to  "a"  or  L/2  for  four-point  or  three-point  beam  systems 
Numerical  factor  (See  Equation  18b) 

Maximum  contact  pressure  at  the  load  application  point 
Radius  of  the  corners  of  the  beam  (See  Figure  5) 

Speed  of  loading 
Time  of  loading 

Vaiiable  beam  distances  (See  Figure  lc) 

Variable  distance  (failure  site  location)  on  either  side  of  the  load 
contact  point  (See  Appendix  D) 

Coordinate  axes  (See  Figure  1) 

Can  also  be  displacement  or  distance  from  neutral  axis  as  defined  in 
the  text 

Beam  curvature  parameters  (See  Equations  3a  and  4a) 

Antic lastic  curvature  factor  (See  Appendix  A) 

A  numerical  factor  associated  with  the  tensile  stress  caused  by  load 
contact  (See  Appendix  D) 

Y  ■  i&(l-v2)/d2p2  (See  Appendix  A) 

Strain  in  the  x,  y,  and  z  directions 

Strain  rate 

Percent  error,  defined  as  [(o.-o  )/ax)100.  A  negative  error  indicates 
the  simple  beam  equations  la  andxlb  underestimate  the  true  stress;  a 
positive  error  is  an  overestimate. 

Angle  of  a  plane  inclined  to  x-axis 

Angle  of  a  plane  inclined  to  the  x-axis  at  which  the  principal  stress 
is  maximum 

Coefficient  of  friction 
Poisson's  ratio 

Radius  of  curvature  of  a  beam  due  to  bending 
Contact  radius  of  a  support  point  (See  Figure  4) 

Contact  radius  of  a  load  point  (See  Figure  4) 

Initial  radius  of  the  curvature  of  a  beam 

Bending  stress  in  a  beam  as  defined  by  simple  beam  theory  or  mean 
fracture  stress 
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Normal  stress  (See  Appendix  C) 

Maximum  principal  stress  (See  Appendix  C) 

Scale  parameter  or  characteristic  value  associated  with  a  Weibull 
analysis 

Stress  in  the  x  direction  (along  the  beam  length) 

Stress  in  the  z  direction  (along  the  beam  width) 

Shear  stress  due  to  torsion  (See  Appendix  C) 

Angle  of  twist  along  the  specimen  length  (See  Figure  3  and  Appendix  C) 
in  Radians 

Angle  of  twist  between  a  pair  of  load  and  contact  points  (See  Figure  3 
and  Appendix  C)  in  Radians 


INTRODUCTION 


There  has  been  an  increase  in  interest  and  activity  in  recent  years  in  both  the 
research  and  development  of  ceramic  materials  ard  their  practical  application  to 
engineering  structures . 

Flexural  testing  is  (and  will  likely  remain)  the  primary  source  of  uniaxial 
strength  data,  either  for  quality  control  or  design  data  purposes.  An  impediment  to 
the  use  of  flexural  strength  data  in  either  application  is  the  lack  of  standard  test 
methods  and  the  presence  of  experimental  error  in  current  practices. 

In  1973,  a  tentative  unapproved  set  of  standards*  was  prepared  by  the  Army 
Materials  and  Mechanics  Research  Center  (AMMRC)  as  it  was  called  at  that  time,  and 
distributed  to  interested  and  involved  organisations.  This  set  of  unofficial  stand* 
ards,  which  included  such  test  methods  as  flexure,  tension,  creep,  stress  rupture, 
fatigue,  and  spin  testing,  was  discussed  at  several  meetings  of  government  and  in* 
dustry  representatives.  A  number  of  worthwhile  suggestions  evolved.  However,  it 
was  apparent  that  these  tentative  standards  were  inadequate  and  thus  not  approved. 
Recently,  however,  interest  was  revived  at  AMMRC,  now  called  the  U.S.  Army  Materials 
Technology  Laboratory  (MTL),  in  finalizing  standard  tests  for  brittle  materials.  It 
was  viewed  that  the  original  tentative  standards,  dated  2  April  1973,  represented 
the  ideal  goal  but  were  far  too  inclusive  to  realistically  establish  testing  require 
ments  which  would  provide  valid  results  at  this  time.  It  was  decided  to  concentrate 
upon  developing  a  standard  method  for  flexural  strength  testing. 

The  objective  of  this  report  is  to  recommend  beam  test  systems  such  that  accu* 
rate  fracture  strength  measurements  will  result  when  testing  brittle  materials 
within  the  elastic  regime.  This  report  differs  from  Reference  1  in  the  following 
ways: 


1.  Discussions  of  a  "Reference  Standard"  beam  system  have  been  deleted  in 
deference  to  a  different  approach  adopted  in  Reference  2. 

2.  The  error  table  for  nonlinear  stress  has  been  eliminated  because  it  was 
redundant  with  errors  analyzed  from  "wedging  stresses." 

3.  The  twisting  error  analysis  has  been  reanalyzed  as  a  plane  stress  condition. 

4.  Additional  analyses,  refinements,  and  corrections  to  the  original  work  have 
been  included  where  appropriate. 

No  attempt  has  been  made  to  determine  the  influence  of  each  error  upon  the 
total  error  of  the  system.  It  is  assumed  that  each  error  is  independent  of  all  the 
other  error  sourcos.  Thus,  for  consistency  and  simplicity,  the  total  error  within 
the  system  is  assumed  to  be  the  sum  of  the  parts.  Errors  in  flexure  testing  of 
beams  are  either  due  to  assumptions  entailed  in  simple  beam  theory,  or  to  sources 
arising  from  external  load  applications.  The  sources  of  error  are  discussed  in  the 
following  sections. 


*  Military  Standards.  Teat  Methods  for  Structural  Ceramics,  2  April  1973. 

1.  BARATTA,  F.  I.  Rtquirtmtnis  /or  Fkxurt  Testing  of  brittle  Material*.  U.  S.  Army  Materials  and  Mechanics  Research  Center,  AMMRC 
TR  82-20,  April  1982,  ADA  1 13937. 

2.  U.  S.  ARMYMILrSTD-1942(MR).  Fhxun  Ttstinf  of  High  Ptrfbrmtmrr  a!  Ambient  T*mp*f*turt.  November  1983. 
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ERRORS  FROM  SIMPLE  BEAM  THEORY  ASSUMPTIONS 


The  rectangular  beaa  configuration  is  attractive  as  a  strength-test  vehicle 
because  of  its  simple  shape  and  apparent  ease  of  load  application*  as  well  as 
analysis  and  reduction  of  data.  Rods  of  circular  cross  section  are  also  used  in 
beam  tests,  but  usually  for  specialised  testing.  Because  a  beam  of  circular  cross 
section  is  not  as  frequently  used  as  the  rectangular  beam,  only  the  rectangular 
cross  section  is  examined  in  the  discussions  to  follow.  Referring  to  Figures  la  and 
2a  for  dimensions  and  applied  loads,  the  simple  beam  formulas  for  maximum  stress  in 
flexure  are: 

■  3P(L-t)/2bd2  for  a  four-point  beam  (la) 

■  3PL/2bd2  for  a  three-point  beam  (lb) 

A  critical  review  of  simple  beam  theory  assumptions  will  yield  ranges  of  geom¬ 
etry  ratios  by  which  the  theory  can  be  validly  applied.  These  assumptions  are 
listed  below,  as  well  as  their  associated  inferences  in  terms  of  an  error  analysis: 

1.  Transverse  planes  perpendicular  to  the  longitudinal  axis  of  the  beam  remain 
plane  after  the  beam  is  deflected. 

2.  The  modulus  of  elasticity  in  tension  is  equal  to  the  modulus  of  elasticity 
in  compression.  Also,  the  beam  matarial  is  isotropic  and  homogeneous. 

3.  The  maximum  deflection  must  be  small  compared  to  the  beam  depth. 

4.  The  beam  must  deflect  normally  under  elastic  bending  stresses  but  not 
through  any  local  collapse  or  twisting. 

5.  Stresses  in  the  longitudinal  direction  are  independent  of  lateral  displace¬ 
ments. 

Each  of  the  above  assumptions  is  examined  in  detail,  where  possible,  so  that  the 
required  rectangular  beam  geometry  ratios  can  be  determined  as  a  function  of  the 
associated  errors. 

Assumptions  1  and  2  together  imply  that  stress  and  strain  are  proportional  to 
the  distance  from  the  neutral  axis,  and  the  stress  does  not  exceed  the  proportional 
limit  of  the  material.  These  assumptions  disregard  the  effect  of  any  shearing 
resistance  and  make  impossible  the  use  of  the  flexure  formula  for  curved  beams  of 
large  curvature. 

Assumption  1  and  the  above  implication  suggest  that  the  bending  stress  is  pro¬ 
portional  to  the  distance  from  the  neutral  axis  to  the  outer  surface  of  the  beam. 

This  assumption  is  valid  if  flexure  of  the  beam  could  be  attained  without  applying 
local  forces  to  the  beam.  However,  practical  flexure  test  systems,  such  as  those 
shown  in  Figures  la  and  2a,  which  utilize  rour-point  and  three-point  beams,  require 
direct  contact  of  the  fixture  to  the  specimen  to  apply  loads  and  thus  moments  to  the 
specimens.  At  the  point  of  contact  there  will  be  compressive  stress  in  the  beam 
depth  direction  resulting  in  a  local  variation  from  linearity  in  the  landing  stress.3 

3.  TIMOSHENKO,  S.,  ud  GOOD1ER,  J.  N.  Theory  off  Jtidty.  2nd  Ed.,  McGnw-MIl  Book  CO.,  Inc.,  Nt»  Yoik,  1951. 
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Figure  1.  Four-point  loading. 


Figure  2.  Threa-point  loading. 
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Because  this  contribution  to  bonding  stress  nonlinearity*  referred  to  as  wedging 
stress*4  is  caused  by  external  load  application*  it  will  be  discussed  further  in 
detail  under  the  section  entitled  Errors  Proa  External  Influences. 

An  error  source  that  is  internal  to  the  beaa  arises  because  of  the  tssuaption 
that  the  aodulus  of  elasticity  in  tension  is  equal  to  that  in  coapression*  By  ■  Bp. 
Chaalis5  has  derived  in  closed  fora  the  solution  for  the  tensile  bending  stress  when 
Ey  i  Ec.  After  soae  aanipulation  of  the  appropriate  formulas,*  the  tensile  stress 
due  to  bending  is  given  by: 

°x  "  (ob/2)l1  +  (VEc)l/2^  (2) 

for  both  the  four-point  and  the  three- point  loaded  beaas.  The  resulting  percent 
error  is  given  in  Table  1. 


Table  1.  ERROR  WHEN  ET  f  Ec 


ET/Ec 

I  Error 

ET/Ec 

*  Error 

0.20 

♦36.2 

1.025 

-0.6 

0.40 

♦22.5 

1.050 

-1.2 

0.60 

♦12.7 

1.075 

•1.6 

0.00 

♦5.6 

1.10 

-2.4 

0.90 

♦2.6 

1.16 

-3.5 

0.925 

♦1.9 

1.20 

-4.6 

0.950 

♦1.3 

1.30 

-6.5 

0.975 

♦0.6 

1.40 

-6.4 

1.00 

0 

1.60 

-11.7 

1.60 

-14.6 

2.0 

-17.2 

Although  the  errors  associated  with  neglecting  to  account  for  anisotropy  and 
nonhoaogeneity  of  the  test  material  are  not  considered  here*  they  are  briefly  men¬ 
tioned  in  the  following  paragraphs  so  that  the  reader  will  be  aware  of  such  possi¬ 
bilities. 

If  the  beam  is  anisotropic,  the  bending  stress  formula  is  exactly  the  same  as 
the  elementary  theory  except  that  the  application  of  a  bending  moment  can  produce 
twisting  moment.  According  to  lekhnitskii ,6  determining  the  accompanying  shear 
stress  produced  by  bending  a  rod  of  rectangular  cross  section*  having  only  one  plane 
of  olastic  symmetry  normal  to  the  axis,  is  very  complicated.  (Composite  and  crystal 
structures  are  excluded  here  as  test  materials.)  If  the  degree  of  anisotropy  for 
ceramic  material  is  slight,  it  may  be  permissible  to  assisse  that  the  error  when 
ignoring  this  effect  on  the  fracture  stress  will  also  be  small. 

Nonhomogeneity  of  the  test  material  infers  variation  of  the  elastic  modulus. 

It  has  been  observed*  that  in  plates  of  hot-pressed  silicon  nitride,  the  modulus  of 
elasticity  at  the  surface  is  several  percent  different  than  that  of  the  center.  This 


♦Private  discussion  with  E.  M.  Lenoe,  MTL.  .  ,  .  ,, 

4.  TIMOSHENKO.  S.  Strength  of  Materials  3rd  Ed..  D.  Van  Most  rand  Co..  Inc.,  N.Y.,  1958,  snd  Part  II.  2nd  Eo.,  D.  Van 
Nostrand  Co.,  Inc.,  N.Y.  1941. 

5.  CHAMUS.  C.  C.  Analysis  of  Three-Point-Beud  Test  for  Materials  with  Unequal  Tension  and  Comprtssiw  Broixriies  NASA  TN  07572, 

6  LEKHNITSKII,  S.  G.  Theory  of  Elasticity  of  an  Anisotropic  Mastic  Body  Hoklcn-Day  Series  in  Mathematical  Physics,  I.  L.  BrandsUtler, 
ed„  1963,  p.  204. 
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is  also  an  area  in  whiyh  further  analysis  will  be  required  to  assess  the  error  appli¬ 
cable  to  four-point  and  three-point  loaded  beams  when  the  modulus  of  elasticity 
varies  through  the  vcam  thickness. 

If  a  rectangular  beam  has  initial  curvature  pc,  the  error  can  be  determined 
from  an  analysis  pro-  ided  by  Timoshenko. 4  The  general  bending  stress  ox  in  a  curved 
beam  due  to  a  pure  moment  is  given  by  the  following: 

°x  *  «cCVbdPc}  C5) 


where 


(d/2pc)  -  (ec/pcl 
c  (sc/pc)[1  -  (d/2pc)] 

ec/pc  "  [fd/Pc)2/l2][1  +  Cd/pc)2/lS]. 


(3a) 

(3b) 


Since  the  bending  stress,  according  to  simple  beam  theory,  is  =  6M^/bd2,  and 
putting  in  the  same  terms  as  (3)  above,  we  have: 


°b  =  “bV^c5’ 


(4) 


where 


“b  *  6(pc/d) .  (4a) 

The  percent  error  F  for  a  beam  of  rectangular  cross  section  and  of  beam-to-depth 
initial  curvature  p  /d  resulting  in  a  neutral  axis  shift  of  e  /P  is: 

C  CO 

e  *  100[(ob  -  ac)/ac].  (5) 

The  resulting  error  for  a  beam  of  rectangular  cross  section  bent  by  a  pure  moment  as 
obtained  from  (5)  is  given  in  Table  2  as  a  function  of  initial  curvature.  It  is  as¬ 
sumed  that  an  analogous  analysis  applied  to  a  three-point  loaded  beam  would  produce 
similar  results. 

The  validity  of  the  assumption  that  the  strain  is  proportional  to  the  distance 
from  the  neutral  axis  and  that  stresses  are  independent  of  lateral  displacements  is 
dependent  upon  the  ratio  of  the  beam  width,  to  its  depth.  Anticlastic  curvature  of 
rectangular  beams  or  plates  with  intermediate  ratios  of  b/d  can  lead  to  erroneous 
results  using  simple  beam  theory;  see  Timoshenko.7  Of  course,  if  the  beam  can  be 
considered  infinite  in  wi.d.h,  like  a  plate,  the  correction  of  the  bending  stress  is 
simply8  1/(1  -  v2).  The  question  arises  as  to  what  ratios  of  b/d  are  appropriate 
for  the  application  of  simple  beam  theory.  Ashwell9  examined  in  detail  the  anti¬ 
clastic  curvature  of  rectangular  beams  and  plates  and  provided  the  answer  to  this 


7.  TIMOSHENKO,  S.  Letter  to  the  Editor.  Mechanical  Engineering,  v.  45,  no.  4,  April  1923,  p.  259-260 

8.  B  ARATTA,F.  I.  When  Is  a  Beam  a  Piute?  J.  Amer,  Cer.  Soc.,  v.  64,  no.  5, 1981,  p.  c-86 

9.  ASHWELL,  D.  G.  The  Anticlastic  Curvature  of  Rectangular  Beams  and  Plates.  J.  Roy.,  Aero.  Soc.,  v.  54, 1950,  p. 
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question.  The  pertinent  formulas  taken  from  Reference  9  are  given  in  Appendix  A. 
These  equations  were  applied  to  ceramic  materials  with  Poisson's  ratio  v  equal  to 
0.25  and  the  ratio  of  Young's  modulus  to  fracture  stress  E/o^  of  1000  to  determine 
the  percent  error*  using  simple  beam  theory  as  a  function  of  b/d  which  is  shown  in 
Table  3. 


Table  2.  ERROR  CAUSED  BY 
INITIAL  BEAM  CURVATURE 


Pc/4 

X  Error 

1 

35.1 

2 

16.7 

3 

10.3 

4 

8.4 

10 

3.2 

15 

2.2 

20 

1.7 

40 

0.8 

100 

0.3 

l  * 

100  [(ab-ac)/<*c3 

Note:  All  errors  are  negative. 


Table  3.  ERROR  CAUSED  BY  EFFECT 
OF  ANTICLASTIC  CURVATURE 


E/ob  -  1x10* 


b/d 

%  Error 

1.0 

0 

15.0 

0 

20.0 

0.1 

30.0 

0.6 

40.0 

1.5 

50.0 

2.6 

100.0 

4.7 

500.0 

5.9 

1000.0 

6.1 

00 

+  (-v*)100  -  -6. 25* 

Note:  All  errors  are  negative. 


if  the  maximum  deflection  is  not  small  compared  to  the  beam  depth,  linear  beam 
theory  cannot  be  employed  without  an  error.  West10  examined  large  deflections  of 
three-point  loaded  beams,  and  from  such  results  a  definitive  ratio  of  beam  length- 
to-depth  can  be  determined  for  valid  application  of  simple  beam  formulas.  Since  for 
most  brittle  materials  values  of  E/o^  range  from  approximately  500  to  1000,  the 
former  value  was  used  to  compute  the  percent  error  because  it  would  yield  the  largest 
error.  Although  the  analysis  was  applied  to  a  three-point  loaded  beam,  the  method 
was  extended  to  determine  errors  for  four-point  loaded  beams  as  well.  The  results 
of  the  calculations  using  the  mentioned  analysis10  are  presented  in  Table  4,  which 
gives  errors  for  four-point  and  three-point  loaded  beams  as  a  function  of  L/d. 


Table  4.  ERROR  FOR  BEAMS  WITH  LARGE 
DEFLECTION 


E/ob  -  500 


L/d 

Four- Point 

Three-Point 

0 

0 

0 

25 

0.1 

0.1 

50 

0.6 

0.4 

100 

1.4 

1.0 

150 

2.5 

1.8 

200 

4.1 

2.9 

250 

7.0 

4.9 

Note:  All  errors  are  negative. 


•Ashwell  considered  a  beam  bent  by  a  constant  moment  analogous  to  the  four-point  beam  loading  cue,  which  should  represent  a 
conservative  bound  on  b/d  for  the  three-point  beam,  u  well. 

10.  WEST,  D.  C.  Flexure  setting  of  Plastics.  Exp.  Mech.,  v.  21,  no.  2,  July  1964. 
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It  is  implicit  in  the  assumptions  given  in  Reference  10  that  the  loads  and 
moments  are  applied  to  the  beam  in  an  ideal  manner  with  no  friction  occurring  between 
the  load  application  points  and  the  beam.  Ritter  and  Wilson11  have  determined  a  beam 
length- to-depth  limit  based  on  the  minimization  of  friction  effects  when  large  de¬ 
flections  occur.  The  friction  effect  considered  is  that  which  gives  rise  to  a  moment 
caused  by  the  slope  at  the  lead  application  point.  Not  considered  in  the  analysis11 
are  the  effects  of  friction  due  to  a  moment  acting  out  of  the  neutral  plane  of  the 
beam,  lateral  contraction  or  extension,  and  changes  in  moment  arms  due  to  contact 
point  tangency  shift.  These  factors  will  be  discussed  later. 

Returning  to  the  results  of  Reference  11,  an  inequality  for  the  four-point- 
loaded  beam  which  provides  a  limit  is  given  in  the  following: 

(L/d  -  a/d)/(E/ob)  <0.3  (6) 

to  insure  negligible  nonlinear  deflections  and  friction  effects.  The  value  of  0.3 
was  obtained  from  limiting  the  slope  to  less  than  15°  between  the  beam  in  the  loaded 
and  unloaded  positions  at  the  outermost  support  point.  If  the  minimum  value  of  E/ob 
is  chosen  to  be  500,  then  we  determine  that  for  a  four-point  loaded  beam  (6)  becomes: 

L/d  -  a/d  <  150.  (7) 

It  is  noted  from  Table  4  that  neglecting  beam  deflections  resulted  in  greater 
error  in  calculation  of  bending  stress  for  the  four-point  loaded  beam  than  for  the 
three-point  loaded  beam.  For  conservatism,  therefore,  it  will  be  assumed  that  (7) 
is  applicable  to  the  three-point  loaded  beam  as  well,  with  a/d  *  0.  Thus  (7)  becomes 

L/d  <  150.  (8) 

It  appears  that  these  limits  are  compatible  with  those  values  given  in  Table  4  such 
that  reasonable  L/d  ratios  can  be  chosen  that  will  result  in  small  errors  when  mini¬ 
mizing  deflection. 

One  of  the  last  requirements,  no  buckling  of  the  beam,  is  easily  fulfilled  for 
ceramic  materials  with  beam  dimensions  of  practical  test  configurations.  The 
reader  can  readily  verify  this  statement  by  referring  to  Timoshenko  and  Gere.12 

Accuracy,  which  is  inferred  in  the  above  restrictions,  is  also  dependent  upon 
the  manner  of  load  application,  beam  geometry,  loading  fixtures,  and  surface  prepara¬ 
tion.  Although  specimen  size  will  not  affect  accuracy  except  for  extremely  small 
geometries,  it  will  alter  the  magnitude  of  the  stress  level  at  failure,  and  this 
must  also  be  considered.  These  subjects  are  discussed  in  the  following  paragraphs, 
^nd  guidelines  for  specimen  geometry  and  minimization  of  errors  are  provided. 

First  to  be  considered,  however,  are  the  merits  of  a  four-point  beam  loading 
system  as  compared  to  the  three-point  beam  loading  system. 


11.  RITTER,  J.  E.,  and  WILSON,  W.  R.  D.  Friction  Effects  in  Four-Point  Bending,  ASLE  Transactions,  v.  18,  no.  2,  p.  130-134,  presented  at 
the  29th  Annual  Meeting,  April  28-May  2, 1974. 

12.  TIMOSHENKO,  S.,  and  GERE,  J.  M.  Theory  of  Elastic  Stability.  McGraw-Hill  Book  Co.,  Inc.,  New  York,  1961. 
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FOUR-POINT  AND  THREE-POINT  LOADING 


The  bending  moment,  from  which  the  desired  fracture  stress  is  computed  in  an 
idealized  four-point  beam  loading  system,  as  shown  in  Figure  la,  is  constant,  and 
there  are  no  horizontal  or  vertical  shear  stresses  within  the  inner  span.  However, 
the  bending  moment  in  an  idealized  three-point  beam  loading  system,  shown  in  Fig¬ 
ure  2a,  is  linearly  dependent  upon  the  distance  from  the  nearest  support  to  the 
fracture  origin,  and  thus  requires  an  additional  distance  measurement  to  determine 
the  fracture  stress.  Also,  the  shear  stresses  for  the  three-point  beam  loading  sys¬ 
tem  are  developed  over  the  full  span,  thus  deviating  from  the  ideally  sought  uni¬ 
axial  stress  state  present  in  the  four-point  beam  loading  system. 

Wedging  stresses4  occur  under  all  points  of  load  application  during  flexure 
testing  of  beams.  The  effect  of  the  wedging  stress  occurring  at  the  inner  load 
points  of  a  four-point  beam  test  is  to  cause  a  deviation  from  the  idealized  calcu¬ 
lated  constant  stress  at  the  two  local  regions.  However,  if  the  ratio  of  half  the 
distance  between  the  outer  span  L  and  inner  span  &,  called  a,  to  beam  depth  d  is 
great  enough,*  the  stress  reduction  will  not  only  be  small  but  will  decay  rapidly, 
and  the  stress  predicted  by  simple  beam  theory  will  be  developed.  Yet,  the  maximum 
stress  computed  by  simple  beam  formula  for  the  three-point  beam  system  is  never 
attained.  The  actual  maximum  stress  occurs  at  a  short  distance  either  side  of  the 
center  of  the  load  application  point,  which  can  cause  fracture  at  these  sites, 
rather  than  at  the  center,  according  to  Rudnick  et  al.13  This  observation  has  also 
been  confirmed  by  Oh  and  Finnie,14  where  only  for  a  material  with  no  scatter  in 
strength  will  the  fracture  location  of  a  three-point  loaded  beam  be  theoretical ly+ 
located  at  the  center  load  point. 

Brittle  materials  are  affected  by  size.  Compensation  can  be  realized  through 
the  use  of  statistical  analysis  offered  by  Weibull.15  Although  the  four-point  beam 
system  assures  a  simple  stress  state  which  is  easier  to  analyze15  than  the  more  com¬ 
plex  biaxial  stress  state  associated  with  the  three-point  beam  specimen,  this  will 
be  less  of  a  consideration  if  the  beam  is  designed  properly.  Nevertheless,  the 
three-point  loaded  beam  system  is  preferred  when  investigating  material  or  process 
development,  because  of  smaller  specimen  size,  or  when  attempting  to  pinpoint 
fracture  origin  location. +  On  the  other  hand,  the  four-point  loaded  beam  is  pre¬ 
ferred  when  determination  of  strength  for  design  purposes  is  desired,  because  the 
center  span  is  uniaxially  stressed,  i.e. ,  no  shear  stresses  exist.  It  is  concluded 
that  each  of  these  systems  is  suited  for  a  particular  application  and  each  has  dif¬ 
ferent  advantages  and  disadvantages. 

Each  of  these  beam  systems  will  be  subjected  to  external  influences  which  will 
affect  the  accuracy  of  the  test  results.  These  external  influences,  directly  or 
indirectly  caused  by  the  application  of  loads  through  the  test  fixtures,  will  lead 
to  either  configuration  constraints  or  errors. 


•This  requirement  will  be  discussed  subsequently. 

tin  Reference  14,  the  authors  considered  only  a  statistical  analysis  and  ignored  wedging  stress  considerations, 
f  Private  communication  with  R.  W.  Rice  of  N.  R.  L. 

13.  RUDNICK,  H„  MARSCHALL,  C.  W„  DUCKWORTH,  W.  H.,  and  ENRICK,  B.  R.  The  Evaluation  and  Interpretation  of  Mechanical  Prop¬ 
erties  of  Brittle  Materials,  AFME  TR  67-316,  April  1968. 

14.  OH,  H.  L,  and  FINNIE,  I.  On  the  Location  of  Fracture  tn  Brittle  Solids  •  l.  Due  to  Static  Loading.  Int.  J.  of  Fracture  Mechanics,  v.  6, 
no.  3,  September  1970,  p.  287-300. 

15.  WEIBULL,  W.  Statistical  Theory  of  Strength  of  Materials.  Royal  Swedish  Institute  for  Engineering  Res.,  Proc.  no.  151, 1939,  p.  1-45. 
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ERRORS  FROM  EXTERNAL  INFLUENCES 


The  maj or  influence  on  the  accurate  determination  of  flexure  strength  of  a  bean 
in  bending  arises  from  the  application  of  load  through  the  fixtures  to  the  specimen. 
The  idealizations  indicated  in  Figures  la  and  2a  are  rarely  met,  and  usually  tests 
are  conducted  using  a  convenient  rigid  loading  head  and  support,  member  as  depicted 
in  Figures  lb  and  2b.  The  constraints  on  either  the  loading  fixture  or  the  specimen 
and/or  errors  resulting  from  such  fixture  designs  are  many.  Such  constraints  or 
errors,  which  are  discussed  in  turn,  are  caused  by: 

1.  eccentric  loading 

2.  span  dimensions 

3.  beam  twisting 

4.  friction 

5.  contact  stresses 

6.  wedging  stresses 

7 .  beam  overhang 

8.  contact  point  tangency  shift 

9.  specimen  preparation 

10.  load  readout 

11.  specimen  dimension  measurement 

Eccentric  Loading 

a.  Four- Point  Loaded  Beams 

When  calculating  bending  stress  by  simple  beam  theory  formula  for  four-point 
loaded  beams,  it  is  usual  to  assume  that  the  moment  within  the  inner  span  l  is  con¬ 
stant.  However,  if  a  loading  head  that  can  only  translate  is  used,  as  idealized  in 
Figure  lb,  it  is  impossible  to  attain  this  idealized  moment  condition  when 
Xi  /  X3  -  x^;13»16  this  is  shown  in  Figure  lc.  The  ratio  of  cx/a^>  from  Appendix  B, 
is: 

ax/°b  =  - - -  *i/a  •  (9) 

L  (Pz  +  Pa)/?  . 

The  loads  and  distances  are  also  shown  in  Figure  lc,  and  a  is  the  value  of  a*  with 
perfect  load  location.  The  error  is  magnified  by  the  ratio  of  xj/a.  (Of  course,  if 
Pi  «  p2  *  P3,  which  implies  exact  location  of  the  points  of  load  application,  there 
is  no  error.)  In  order  to  estimate  the  magnitude  of  such  an  error  it  was  assumed  in 

16.  HOAGLAND,  R.  G.,  MARSCHALL,  C.  W.,  and  DUCKWORTH,  W.  H.  Reduction  of  Errors  in  Ceramic  Bend  Tests.  J.  Amer.  Cer.  Soc., 
y.  59,  no.  5*6,  May-June  1976,  p.  189*192. 
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Appendix  B  that  the  upper  two  load  points  in  Figure  lc  were  at  a  fixed  distance 
x2  “  xl  *  *■  and  were  constrained  to  translate  vertically  during  loading,  and  that 
the  loading  head  would  be  located  such  that  Xi  +  X3  -'  x2.  This  method  of  loading, 
being  the  most  convenient,  is  often  adopted  by  many  investigators,  and  therefore  the 
resulting  error  determination  is  not  unrealistic. 

The  analysis  was  accomplished  by  simply  enforcing  the  condition  that  the  dis¬ 
placement  at  xi  must  be  equal  to  the  displacement  at  x2  in  the  deflection  equation. 
This  results  in  the  following  relationships  between  ox  and  0,  in  terms  of  the  load 
eccentricity  ratio  e/L  (see  Appendix  B  for  details): 


[(e/L  +  a/L) / (a/L) ] [ 1- (e/L  +  a/L)  -  A/L] { (A/L) [2- (e/L  +  a/L)]-2[l-(e/L  +  a/L)]2} 

ax/0  a  - - 

3 (e/L  +  a/L) [1-  A/L  -(e/L+  a/L)]-(l-  A/L)2 

(10) 

where  the  parameter  a  defined  as  (L-A)/2,  e  defined  as  (ai-a);  and  A  and  L  are  shown 
in  Figure  1. 

Most  workers  in  the  testing  field  utilize  either  a  1/3-point  (a/L  =  1/3  and 
A/L  =  1/3)  or  a  1/4-point  (a/L  =  1/4  and  A/L  =  1/2)  loading.  Thus  by  substitution  of 
these  parameters  into  Equation  10,  we  obtain: 


(ox/ob) 


[3(e/L)+l] (1/3  -  e/L) [ 1/3 (5/3  -  e/L)-2(2/3  -  e/L)2] 
[3(e/L)+l] (1/3  -  e/L)  -  4/9 


(°x/ab^/L=]  - 


[4 (e/L) +1] (1/4  -  e/L) [1/ 2(7/4  -  e/L)-2(5/4  -  e/L)2] 
[3 (e/L)  +  3/4] (1/4  -  e/L)  -  1/4 


(ID 


(12) 


The  reader  is  cautioned  that  for  given  values  of  A/L  there  exists  a  limit  on 

e/L  in  (10),  (11),  and  (12);  that  is,  if  aj  is  such  that  either  P2  or  P3  =  0,  the 

test  system  changes  from  four-point  to  an  eccentric,  three-point  loading  (see 
Appendix  B) ,  and  the  above  equations  become  invalid. 

The  error,  defined  as  [(ob  -  cO/OjJlOO,  was  determined  from  (11)  and  (12)  for 
the  1/3-point  and  1/4-point  loaded  beams  and  is  shown  in  Tables  5  and  6  as  a  function 
of  e/L.  Only  negative  values  of  e/L  were  considered  in  (11)  and  (12)  because  when 
e/L  <  0,  a*  >  ob.  A  negative  value  of  e/L  corresponds  to  the  inner  load  bearing 
which  is  offset  closer  to  the  outer  load  bearing  (aj  <  a) .  An  error  of  similar 
magnitude,  but  larger  and  of  opposite  sign,  exists  at  the  other  inner  load  bearing, 

which  is  why  Tables  5  and  6  show  +_ values.  Tables  5  and  6  show  that  for  correspond¬ 

ing  e/L,  when  aj/L  /  a2/L,  the  1/3-point  loading  system  results  in  lesser  error  than 
the  1/4-point  loading  system.  Also,  in  accordance  with  the  above  discussion,  e/L 
in  Tables  5  and  6  is  limited  to  a  range  of  +0.0443  and  +0.0465.  The  errors  indicated 
in  these  tables  can  be  minimized  by  designing  the  loading  fixture  so  that  the  inner 
and  outer  spans  are  independently  fixed.  Also.,  the  inner  span  should  be  designed 
with  accurate  location  adjustment  and  allowed  to  pivot  as  recommended  by  Hoagland 
et  al.16 
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Table  5.  ERROR  DUE  TO  ECCENTRIC  LOAD  APPLICATION 
FOR  A  1/3-F0UR-P0INT  LOADED  BEAN.  FOR  EITHER  A 
NON-PIVOTING  OR  PIVOTING  LOADING  HEAD 


Whan  t/L  -  1/3  and  at/L  /  a2/L 


e/L  »  +(«i/L  -  1/3) 

NON-PIVOTING 
♦  X  Error 

PIVOTING 
♦  X  Error 

0.0 

0.0 

0.0 

0.0010 

0.7 

0.1 

0.0019 

1.3 

0.2 

0.0038 

2.6 

0.4 

0.0057 

3.8 

0.6 

0.0076 

4.9 

0.7 

0.0095 

6.0 

0.9 

0.0114 

7.0 

1.1 

0.0133 

8.1 

1.2 

0.0333 

16.1 

2.6 

0.0433 

18.7 

3.1 

0.0443 

18.9 

3.2 

Table  6.  ERROR  DUE  JO  ECCENTRIC  LOAD  APPLICATION  . 
FOR  A  1 /4-FQUR” POINT  LOADED,  BEANt  FOR  EITHER  A 
NON-PIVOTING  OR  A  PIVOTING  LOADING  HEAD 

When  t/L  ■  1/2: and  at/L  A  aj/L 


e/L  -  +{ai/L  -  1/4) 

NON-PIVOTING 
+  1  Error 

pivoting 

+  X  Error 

0.0. 

0.0 

0.0 

0.0010 

1.0 

0.2 

0.0020 

2.1 

0.4 

0.0040 

3.8 

0.8 

0.0080 

7.1 

1.5 

0.0120 

1Q.0 

2.2 

0.0160 

12.6 

2.9 

0.0200 

14.7 

3.6 

0.0240 

16.6 

4.2 

0.0280 

18.3 

4.7 

0.0320 

19.8 

5.3 

0.0340 

20.8 

5.6 

0.0400 

21.8 

6.3 

0.0465 

22.9 

7.0 

Many  flexure  fixtures  do  permit  the  loading  head  to  translate  and  pivot.  The 
eccentric  loading  error  in  this  instance  will  be  due  to  the  actual  moment  being 
different  from  the  assumed  moment.  If  the  beam  deflection  is  small,  the  angular 
rotation  of  the  loading  head  can  be  ignored  and  the  maximum  bending  stress  can  be 
determined  utilizing  force  and  moment  equilibrium: 


o 


X 


a 


b 


2e  <e  2e2 
L  a  ~  aL 


(13) 


The  maximum  stress  will  exist  tinder  the  inner  load  bearing  which  is  offset  to 
give  a  larger  moment  arm  (a).  A  similar  error  (for  e/L  <,  0.01),  but  larger  and  of 
opposite  sign,  will  exist  at  the  opposite  inner  load  bearing.  Substituting  for 
either  1/3  or  1/4-point  loading: 


6el 

L2 


8ei 

L2 


(14a) 


(14b) 


The  latter  expression  was  also  derived  by  Jayatilaka. 17  These  errors  at  the 
point  of  maximum  stress  are  also  given  in  Tables  5  and  6  for  comparison.  It  is 
evident  that  a  translating  and  pivoting  loading  head  is  preferred  to  a  rigid  loading 
head  because  the  errors  are  appreciably  less.  This  finding  is  consistent  with 
recommendations  of  Hoagland  et  al.16 


17.  JAYATILAKA,  A.  DeS.  Fnctun  of  Engineering  Brittlt  Mattrlah.  Applied  Science  Publiihen,  LTD,  London,  1979,  p.  187. 
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b.  Three- Point  Loaded  Beans 


The  ratio  of  o„  to  o.  is: 

x  o 


The  percent  error  as  a  function  of  e/L  is  given  in  Table  7.  Notice  that  the  percent 
errors  in  Table  7  are  always  positive,  and  when  the  load  application  point  is  mis¬ 
placed,  such  errors  are  much  less  than  those  of  equivalent  e/L  values  shown  in 
Tables  5  and  6  for  the  four-point  loaded  beams.  Notice  also  that  when  ^e/L  ■  0.500, 
the  error  is  infinite,  i.e.,  the  three-point  loading  model  is  no  longer  valid. 


Table  7.  ERROR  DUE  TO  ECCENTRIC  LOAD 
APPLICATION  FOR  A  THREE-POINT 
LOADED  BEAM 


When  ti/L  /  a2/L  i  1/2 
e/L  -  1/2  -  a^L 


ie/L 

X  Error 

0 

0 

0.025 

0.25 

0.050 

1.0 

0.075 

2.3 

0.100 

4.2 

0.150 

9.9 

0.200 

19.0 

0.250 

33.3 

0.300 

56.3 

0.400 

177.8 

0.450 

426.3 

0.500 

•• 

Note:  All  errors  are  positive. 

Span  Dimensions 

a.  Four- Point  Loaded  Beams 

An  additional  mislocation  error  may  exist  if  the  inner  bearing  span  (£)  or  the 
outer  bearing  span  (L)  are  not  their  prescribed  values,  even  if  they  are  properly 
centered  with  respect  to  each  other.  This  will  alter  the  moment  arm  (a) .  Assuming 
the  inner  span  is  actually  l  *  eg  and  the  outer  span  is  L  -  eg,  then  the  ratio  of 
(j  to  a.  is: 

X  D 


5T’  1  -TVa  '  1)1 

D 


(16) 


where  e_  is  the  error  of  the  inner  and  outer  span  dimensions.  A  similar  error 
(for  e/L  <_  0.01),  but  of  opposite  sign  exists  if  the  inner  span  is  l  -  e,  and  the 
outer  span  is  L  +  e  .  Errors  are  tabulated  in  Table  8  for  the  1/3  and  1/4-point 
loaded  beams.  The  largest  error  magnitude,  occurring  when  the  outer  span  is 
L  -  e  ,  is  reported  in  these  tables. 
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TabH  8.  ERROR  DUE  TO  MRQN8  SRMIS 


±*,/L 

1/3-Four-Point 

♦  *  Error 

1/4-Four-Rolnt 

Three-Point 

0 

0 

0 

0 

0.001 

0.3 

0.4 

0.1 

0.002 

0.6 

0.8 

0.2 

0.005 

1.5 

2.0 

0.5 

0.010 

.3.1 

4.2 

1.0 

0.015 

4.7 

6.4 

1.5 

0.020 

6.4 

8.7 

2.0 

0.025 

8.1 

11.1 

2<6 

0.050 

17.6 

25.0 

5.3  . 

Nota:  All  errors  ara  either  positive  or  nefttlve. 


b.  Three^ Point  Loaded  Beans 


A  simple  analysis  shows  that  if  the  support  span  is  actually  L  -  e  ,  then: 

s 


(17) 


The  error  in  detexmining  the  stress  is  given  in  Table  8.  If  the  support  span 
is  L+e  ,  a  similar  error  occurs  but  it  is  slightly  less  and  of  opposite  sign.  A 
comparison  of  Equations  16  and  17  shows  that  the  four-point  configuration  amplifies 
the  span  error,  whereas  the  error  in  computing  the  stress  for  a  three-point  beam  is 
nearly  the  same  as  the  sp^n  error. 

Beam  Twisting 

A  net  torque  can  result  from  line  loads  being  nonunifoxm  or  nonparallel  between 
pairs  of  load  contact  points  or  if  the  cross  section  of  the  specimen  is  skewed  over 
its  length.13*16 

Such  a  skewed  condition  is  shown  schematically  in  Figure  3  for  a  four-point 
bending  specimen.  The  error  due  to  twisting  has  been  estimated  for  plane  strain 
and  plane  stress  conditions  by  examining  the  maximum  principal  stress  due  to  banding 
and  torsion  and  comparing  it  to  the  bending  stress.16  "Bottoming"  of  the  specimen 
on  the  fixture  was  not  considered.  Bottoming  occurs  when  the  bearing  rollers  contact 
the  specimen  across  its  full  width.  For  the  sake  of  completeness,  bottoming  is  con¬ 
sidered  in  the  analysis  given  in  Appendix  C.  The  maximum  principal  stress,  assuming 
a  plane  strain  condition,  is  derived  in  Reference  1.  The  plane  strain  criterion 
leads  to  slightly  higher  error  estimates,  but  the  plane  stress  criterion  is  more 
appropriate.  The  plain  stress  solution  is  also  given  in  Reference  16*  but  the 
analysis  has  been  extended  in  this  report  to  incorporate  the  case  where  the  spiecimen 
bottoms  on  the  fixture. 

The  maximum  principal  stress  for  either  a  skewed  four-point  or  three-point  beam 
in  bending,  considering  a  plane  stress  condition,  is  given  by: 

0  -  ob/2(l  +  (l/3k2)[(nb/t»)2  ♦  9k22])l/2  (18a) 
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where  ob  is  the  apparent  bend  strength  end  V  is  either  equal  to  "a"  for  four-point 
bending  or  equal  to  L/2  for  three-point  bending.  Also: 

n  -  [3ki(e/ob)/(l  ♦  v)][(d/Ur)^s  +  (d/t'Hp]  (-£-)  (18b) 

where  for  Case  I:  n  ■  1,  failure  occurs  prior  to  bottoming  of  the  specimen  in  the 
loading  fixture,  and  for  Case  II:  n  <  1,  failure  occurs  after  bottoming. 

tt 


(b)  End  View  of  Specimen. 
The  applied  loads  contact 
the  specimen  edges  due  to 
the  twist. 


(c)  End  View  of  Fixture  Showing  +F,  the  Fixture  Twist 
Angle  Between  a  Pair  of  Contact  and  Load  Bearings 


P  P 


Figure  3.  Twitting  of  •  four-point  beam  specimen. 

The  factors  kj  and  k2*  obtained  from  Reference  3  and  given  in  Table  9,  are 
numerical  values  associated  with  the  torsional  stress  component  which  are  dependent 
on  the  ratio  of  b/d.  The  measured  angle  of  twist  (or  skew  angle)  along  the  total 
length  Lt  of  thi  specimen  is  +.  (see  Figure  3c),  and  along  the  Tixture  from  one  sup¬ 
port  point  to  the  adjacent  load  point  is  $p. 

The  mayfmniB  principal  stress  as  given  by  (18a)  can  be  utilized  to  determine  the 
percent  error  for  various  ratios  of  n,  t'/b,  and  b/d.  This  was  accomplished  and  is 

shown  in  Table  10.  Notice  that  the  range  of  n  varies  from  0.20  to  1.00.  It  is 

expected  that  if  bottoming  dous  not  occur  prior  to  fracture  because  of  an  excessive 

twist  angle,  the  maximum  ratio  of  n  that  can  be  attained  is  1.0  and  thus  the  tables 

do  not  accommodate  n  >  1.0. 
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Tabic  9.  k,  AND  kt 


b/d 

*i 

1.0 

0.1406 

0.208 

1.2 

0.166 

0.219 

1.5 

0.196 

0.231 

2.0 

0  229 

0.246 

2.5 

0.249 

0.256 

5.0 

0.291 

0.291 

10.0 

0.312 

0.312 

•b 

0.333 

0.333 

Table  10.  ERROR  OUE  TO  BEAN  TWISTING  PLANE  STRESS  ASSUMPTION* 


v  -  0.25 


t'/b 

b/d 

1.0 

1.2 

1.5 

2.0 

2.5 

5.0 

n  ■  0.20 

1.0 

2.44 

2.22 

2.00 

1.77 

1.62 

1.27 

2.0 

0.63 

0.57 

0.52 

0.46 

0.42 

0.32 

2.5 

C.41 

0.37 

0.33 

0.29 

0.27 

0.21 

5.0 

0.10 

0.09 

0.08 

0.07 

0.07 

0.05 

10.0 

0.03 

0.02 

0.02 

0.02 

0.02 

0.01 

•• 

0 

0 

0 

0 

0 

0 

n  ■  0.40 

1.0 

8.57 

7.87 

7.18 

6.44 

5.93 

4.75 

2.0 

2.44 

2.22 

2.00 

1.77 

1.62 

1.27 

2.5 

1.59 

1.44 

1.30 

1.15 

1.05 

0.82 

5.0 

0.41 

0.37 

0.33 

0.29 

0.27 

0.21 

10.0 

0.10 

0.09 

0.08 

0.07 

0.07 

0.05 

«• 

0 

9 

0 

0 

0 

0 

n  ■  0.60 

1.0 

16.20 

15.05 

13.09 

12.63 

11.74 

9.61 

2.0 

5.19 

4.73 

4.29 

3.82 

3.51 

2.78 

2.5 

3.44 

3.13 

2.83 

2.52 

2.30 

1.82 

5.0 

0.91 

0.82 

0.74 

0.65 

0.60 

0.47 

10.0 

0.23 

0.21 

0.19 

0.16 

0.15 

0.12 

m 

0 

0 

0 

0 

0 

0 

n  ■  0.80 

1.0 

23.81 

22.36 

20.87 

19.20 

18.01 

15.09 

2.0 

8.57 

7.87 

7.18 

6.44 

5.93 

4.75 

2.5 

5.82 

5.32 

4.83 

4.31 

3.95 

3.14 

5.0 

1.59 

1.44 

1.30 

1.15 

1.05 

0.82 

10.0 

0.41 

0.37 

0.33 

0.29 

0.27 

0.21 

m 

0 

0 

0 

0 

0 

0 

n  •  1.0 

1.0 

30.74 

29.12 

27.43 

25.50 

24.11 

20.61 

2.0 

12.32 

11.38 

10.44 

9.42 

8.72 

7.06 

2.5 

8.57 

7.87 

7.18 

6.44 

5.93 

4.75 

$.0 

2.44 

2.22 

2.00 

1.77 

1.62 

1.27 

10.0 

0.63 

0.57 

0.52 

0.46 

0.42 

0.32 

«• 

0 

0 

0 

0 

0 

0 

Note:  All  errors  are  negative. 

*An  error  table  based  upon  plain  strain  conditions  Is  In  Reference  1. 
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Friction 


It  has  already  been  shown  in  Table  4,  for  the  two  beast  systems  considered,  that 
the  erro*~  due  to  deflection  will  be  negligible  if  L/d  5  25,  It  appears  that  this 
limit  is  »*ell  within  an  attainable  realistic  geometry  ratio.  Therefore,  the  fric¬ 
tion  effect  at  the  load  and  support  points  will  be  minimized  with  respect  to  large 
deflections.  This  also  implies  that  there  will  be  no  effect  from  friction  on  the 
contact  tangency  shift.  (These  factors  will  be  discussed  subsequently.)  However, 
friction  will  cause  a  moment  acting  out  of  the  plane  of  the  beam  that  can  not  be 
ignored.  This  factor  is  considered  in  the  following. 


When  determining  bend  strength  by  simple  beam  theory,  it  is  usual  to  assume  that 
the  supports  and  load  points  are  frictionless,  whereas  in  fact  they  are  not.  The 
presence  of  friction  in  flexure  tests  with  fixed  load  and  support  points  gives  rise 
to  couples  at  such  locations  as  well  as  axial  forces  at  the  neutral  axis  of  the  beam. 
The  net  axial  force  is  relatively  small  and  therefore  is  ignored  here.  However,  if 
the  moment  is  not  corrected  to  account  for  the  couple  in  the  determination  of  flexure 
stress,  an  error  will  result.  Error  equations  adapted  from  the  results*  available 
in  the  literature16**19  are  given  below  for  the  four-point  and  three-point  loading 
systems : 


e  ■ 


and 

T  ■ 100  (17^7) 


(19) 


(20) 


Such  errors  as  defined  by  the  above  equations  can  be  significant,  according  to 
References  16,  19  and  20.  Newnham19  and  Weil20  reported  that  the  experimental  dif¬ 
ference  in  failure  stress  using  rigid  knife  edges  as  compared  to  roller- type  contact 
points  was  as  high  as  12%  for  silicon  nitride  and  13%  for  graphite. 


Contact  Stresses 


Loads  on  bend  specimens  applied  through  knife  edges  or  small-diameter  rollers 
result  in  high  stresses  under  these  line  loads.  High  compressive  contact  stresses 
can  result  and  cause  local  crushing.  (Also,  shear  stress  near  the  locality  of  the 

load  point  can  be  several  times  higher  than  that  predicted  by  beam  theory.) 

Reference  4  gives  equations  for  determining  the  contact  pressure  between  a 
cylinder  (or  roller)  and  a  flat  surface  (see  Figure  4)  as  a  function  of  the  applied 
load,  modulus  of  each  material,  and  the  roller  radius.  If  it  can  be  assumed  that 
the  two  materials  are  identical  and  that  the  allowable  boaring  pressure  or  contact 
pressure  can  be  as  high  as  twice  the  bend  strength  of  the  material,  then  limits  on 
the  roller  radius  for  both  loading  systems  will  result.  For  example,  from  Reference 
4  we  have: 

"Beam  width  constraint  occurs  also  because  of  friction  tniuwns  to  the  beam's  long  axis.  However,  this  effect  (see  Newnham19)  is 
small  and  thus  not  considered  here. 

18.  DUCKWORTH,  W.  H.,  et  al.  Mechanket-Property  Tests  on  Ctnunic  Bodies.  WADC  TR  52-67,  March  1952,  p.  67-70. 

19.  NEWNHAM,  R.C.  Strength  Tests  for  Brittle  Metertoh.  Ftoc.  of  the  British  Car.  Soc.,  no.  25,  May  1975,  p.  281-293. 

20.  WEIL,  N.  A.  Studies  of  Brit  He  Beherknir  of  Ceremic  Meteriets.  ASD  TR  61-628,  Part  II,  April  1962,  p.  38-42. 
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p, 


0.59*te/2bpi 


(21) 


where  p  Is  the  maximum  contact  pressure.  (Note  that  the  roller  radius  can  ba 
either  pi  or  p*. )  However,  we  shall  assume  that  p  <  2ov.  Also  for  four-point 
loading,  oh  ■  6Pa/bd2,  and  for  three-point  loading,  #b  -  (3/2)PL/bd2.  Substituting 
of  into  (17)  and  solving  for  pi/d  we  obtain 

Pi/d  >  7.25  d/a  for  the  four-point  loaded  beam,  and  (22a) 

Pl/d  >_  29.0  d/L  for  the  three-point  loaded  beaa,  (22b) 

where  it  was  assumed  that  E/o.  ■  1000.  Of  course,  if  the  specimen  and  bearing  are 
made  of  different  materials,  and  if  E/d  is  not  1000,  then  further  calculations  are 
required  to  ensure  that  the  ceramic  doesn’t  locally  crush  or  fracture,  or  that  the 
bearing  does  not  permanently  flatten. 


Wedging  Stresses 

Localized  contact  at  the  load  bearings  can  cause  a  more  subtle  problem,  which 
is  referred  to  as  wedging  stresses.  The  effect  of  the  wedging  stress  is  to  provide 
a  substantial  tensile  Stress  contribution  at  the  compressive  side  of  the  beam  adja¬ 
cent  to  the  load  points.  A  net  tensile  stress  can  not  be  created  if  d/2t'  <  1, 
according  to  Reference  16.  More  importantly,  a  tensile  stress  is  added  to  that 
already  present  due  to  beam  bending  at  the  tensile  side  of  the  beam,  thereby  causing 
a  deviation  from  the  assumed  stress  calculated  by  simple  beam  theory. 
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This  problem  is  generally  treated  in  Reference  3  and  particular  results  from 
von  tfrmdh  and  Seewald21  for  a  sinilar  situation  are  used  to  estimate  this  error. 

An  analysis  for  this  error  is  given  in  Appendix  D.  The  resulting  error  determina¬ 
tions  for  four-  and  three-point  loaded  beams  are  given  in  Table  11.  In  the  calcula 
tion  o i  the  errors,  which  are  a  function  of  a/d  or  L/d,  as  well  as  x’/d,  the  com¬ 
puted  corresponds  to  the  failure  site  location  (x’/d) . 


Tibi#  11.  ERROR  DUE  TO  UE06IN6 


*'/d# 

Loading 

0 

0.125 

0.25 

0.375 

0.50 

0.75 

1.0 

1.50 

a.  Four-point. 

mid 

ft 

♦4.7 

•0.5 

•2.8 

-2.1 

-1.4 

-0.7 

-0.3 

t 

1.5 

♦3.1 

-0.3 

-1.9 

-1.4 

-0.9 

-0.5 

-0.2 

0 

2.0 

+2.3 

•0.2 

-1.4 

-1.1 

-0.7 

•0.4 

-0.2 

3 

3.0 

♦1.5 

-0.2 

-1.0 

-0.7 

-0.5 

-0.2 

-0.1 

0 

4.0 

♦1.1 

-0.1 

-0.7 

-0.5 

•0.3 

-0.2 

-0.1 

0 

5.0 

♦0.9 

-0.1 

•0.6 

•0.4 

-0.3 

-0.1 

0 

0 

6.0 

♦0.8 

-0.1 

•0.5 

•0.4 

-0.2 

-0.1 

0 

0 

8.0 

♦0.6 

0 

-0.4 

•0.3 

•0.2 

-0.1 

0 

0 

10.0 

♦0.4 

0 

-0.3 

•0.2 

-0.1 

-0.1 

0 

0 

15.0 

♦0.3 

0 

-0.2 

-0.1 

-0.1 

0 

0 

0 

20.0 

♦0.2 

0 

-0.1 

-0.1 

-0.1 

0 

0 

0 

40.0 

♦0.1 

0 

-0.1 

-0.1 

0 

0 

0 

0 

60.0 

+0.1 

0 

0 

0 

0 

0 

0 

0 

ea 

0 

0 

0 

0 

0 

0 

0 

0 

b.  Thr##-po1nt 

i  id 

ft 

♦21.6 

-2.4 

•18.8 

-25.4 

t 

+ 

t 

t 

1.5 

♦13.4 

-1.4 

-10.4 

-10.2 

-10.1 

t 

t 

t 

2.0 

♦9.7 

-1.0 

-7.2 

•6.4 

-5.3 

-5.5 

T 

t 

3.0 

+6.3 

-0.7 

-4.4 

-3.7 

-2.7 

-1.9 

-1.3 

t 

4.0 

♦4.7 

-0.5 

-3.2 

-2.6 

-1.8 

-1.2 

-0.6 

-0.1 

5.0 

♦3.7 

-0.4 

-2.5 

-2.0 

-1.4 

-0.8 

-0.4 

0 

6.0 

♦3.1 

-0.3 

-2.1 

-1.6 

-1.1 

-0.6 

-0.3 

0 

8.0 

♦2.3 

-0.2 

-1.5 

-1.2 

-0.8 

-0.4 

-0.2 

0 

10.0 

♦1.8 

-0.2 

-1.2 

-0.9 

-0.6 

-0.3 

-0.2 

0 

15.0 

+1.2 

-0.1 

-0.8 

-0.6 

-0.4 

-0.2 

-0.1 

0 

20.0 

♦0.9 

-0.1 

•0.6 

-0.4 

-0.3 

-0.2 

-0.1 

0 

40.0 

+0.4 

0 

-0.3 

•0.2 

-0.1 

-0.1 

0 

0 

60.0 

+0.3 

0 

-0.2 

-0.1 

-0.1 

-0.1 

0 

0 

aa 

0 

0 

0 

0 

0 

0 

0 

0 

V  Is  the  distance  on  #1th#r  sltf#  of  th«  load  contact  point  whtro  fallur#  occurs. 
^Location  Is  at  or  beyond  out#r  span  limit. 


Beam  Overhang 

The  overhangs  of  the  beam  must  be  great  enough  so  that  the  local  stresses  at 
the  beam  support  points  are  not  amplified  due  to  beam-end  effects.  These  stresses 
are  dampened  out  within  a  distance  equal  to  one  beam  depth.21  Thus,  by  allowing 

1^,  >  L  +  2d  (23) 

beam-end  effects  are  avoided. 

21.  VON  KARMAN,  T.,  and  SEEWALD,  F.  Alhandl  Aarodynam,  but.  Tech.  Hochacfauta,  1946,  p.  256. 
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Contact  Point  Tangency  Shift 

Significant  changes  in  span  length  can  occur  in  both  four-point  and  three-point 
loading  systems  if  contact  radii  of  support  and  load  points  are  large  compared  to 
beam  depth.  The  shift  in  point  of  tangency,  as  shown  by  hi  and  h2  in  Figure  4,  is  a 
function  of  the  contact  radii,  specimen  thickness,  and  the  ratio  of  the  modulus  of 
elasticity  to  the  bend  strength.  For  materials  that  behave  elastically,  such  as 
those  considered  here,  the  change  in  tangency  point  and  thus  the  error  arising  be¬ 
cause  of  the  change  in  moment  arm  from  the  ideal  can  be  predicted  mathematically 
for  linear  systems.  This  is  accomplished  and  is  presented  in  Appendix  E.  The 
approach  was  patterned  after  Westwater22  who  corrected  for  span  shortening  but 
ignored  friction  at  the  support  points  of  a  three-point  loaded  beam.* 

In  Appendix  E  the  formulas  are  derived  for  a  four-point  loaded  beam  and  then 
reduced  to  the  special  case  of  a  three-point  loaded  beam.  These  results  are  put  in 
terms  of  error  functions  assuming  the  simple  beam  theory  is  applied  without  correct¬ 
ing  for  span  shortening,  as  in  the  case  of  the  lower  support,  and  span  lengthening 
between  the  upper  loading  points  shown  in  Figure  4. 

The  errors  are  determined  for  four-point  loaded  beams  of  1/3  and  1/4  loading 
points  as  a  function  of  Pi/d  and  pp/d,  and  the  three-point  loaded  beam  as  a  function 
of  Pi/d  only.  These  errors  are  given  in  Table  12,  where  it  was  assumed  that 
E/a  -  1000. 


Table  12.  X  ERROR  DUE  TO  TANGENCY  POINT  SHIFT 
E/ob  «  1  X  10* 


P*/d 


Loading 

pi/d 

1.0 

2.0 

5.0 

10.0 

a.  Four-point, 

1.0 

0.3 

0.4 

0.7 

1.2 

a/L  =■  1/3 

2.0 

0.5 

0.6 

0.9 

1.4 

4.1 

0.9 

1.0 

1.3 

1.8 

6.1 

1.3 

1.4 

1.7 

2.3 

8.2 

1.7 

1.8 

2.1 

2.7 

10.3 

2.1 

2.2 

2.6 

3.1 

b.  Four-point, 

0.67 

0.4 

0.6 

”  1.2 

2.2 

a/L  -  1/4 

1.35 

0.6 

0.8 

1.4 

2.5 

2.7 

1.0 

1.2 

1.8 

2.9 

4.1 

1.4 

1.6 

2.2 

3.3 

5.5 

1.8 

2.0 

2.6 

3.7 

6.9 

2.2 

2.4 

3.1 

4.1 

c.  Three-point, 

1.0 

0.1 

a/L  =  1/2 

2.0 

0.2 

4.0 

6.0 

0.4  1 
0.6  | 

Regardless  of  pj/d  value 

8.0 

0.8 

10.0 

1.0 

Note:  All  errors  are  positive. 


*Westwatar  also  determined  an  approximate  relationship  for  the  horizontal  load  arising  because  of  tangency  shift.  However,  for 
beams  of  small  deflection,  the  error  is  negligible. 

22.  WESTWATER,  J.  W.  Flexure  Testing  of  Plastic  Materials,  hoc.  ASTM,  v.  49, 1949. 
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Specimen  Preparation 

The  flexure  strength  of  each  brittle  material  is  not  only  supersensitive  to  the 
final  surface  finish  because  the  maximum  tensile  stress  occurs  at  the  beam  surface, 
but  is  also  highly  sensitive  to  prior  finish  history.  For  this  reason  it  is  impos¬ 
sible  to  specify  an  optimum  surface  finish  procedure  for  all  brittle  materials,  so 
that  failure  will  be  due  to  inherent  flaws  related  to  the  material  or  material  proc¬ 
essing,  rather  than  an  imposed  defect  resulting  from  the  finish  process.  Indeed,  the 
designer  or  materials  developer  may  not  be  able  to  specify  a  particular  finish  pro¬ 
cedure.  Therefore,  rather  than  attempt  to  dictate  surface  finish  requirements,  it  is 
suggested  that  each  set  of  reported  test  data  results  be  accompanied  by  surface 
finish  history  and/or  material  process  history,  whichever  is  applicable. 

There  are,  however,  several  specific  recommendations  related  to  surface  finish¬ 
ing  procedures  that  can  be  presented.  Corner  flaws  resulting  from  chipping  or  crack¬ 
ing  during  the  grinding  operation  are  sources  of  low-strength  failure.  Rounding  or 
beveling  of  the  comer  as  depicted  in  Figure  5  appears  to  reduce  premature  failure.23 
Since  a  chamfer  will  double  the  number  of  edges,  thus  doubling  the  source  of  flaw 
locations,  rounding  is  preferred.24  Also,  it  is  important  to  grind  the  edges  and 
flat  surfaces24  by  a  motion  parallel  to,  rather  than  perpendicular  to,  the  specimen 
length.  It  is  further  indicated23  that  finishing  of  the  corner  should  be  comparable 
in  all  aspects  to  that  applied  to  the  beam  surfaces. 


Rounded  Corners  Chamfered  Corners 


Figure  5.  Beam  cross  section. 


If  the  corner  radii  or  chamfer  is  small,  the  error  in  ignoring  the  change  in 
moment  of  inertia  will  be  negligible.  The  limiting  ratio  of  corner  radii  or  45° 
chamfer  dimension  to  beam  depth  can  be  determined  from  the  error  analysis  due  to 
neglecting  the  change  in  moment  of  inertia  given  in  Appendix  F.  This  error  in 
determining  flexure  stress,  when  neglecting  comer  radii  or  45°  chamfer,  is  given 
in  Table  13. 


23.  RICE,  R.  W.  Machining  of  Ceramics.  Proc.  of  the  Second  Army  Materials  Technolony  Conference  -  Ceramics  for  High  Performance 
Applications,  I.  J.  Burke,  A.  E.  Gorum,  and  R.  N.  Katz,  ed.,  Brook  Hill  Publishing  Company,  Chestnut  Hill,  Massachusetts,  1974. 

24.  RICE,  R,  W.  The  Effect  of  Grinding  Direction  on  the  Strength  of  Ceramics.  The  Science  of  Ceramics  Machining  and  Surface  Finishing, 

S.  J.  Scheider  and  R.  W.  Rice,  ed.,  NFS  Special  Publication  348,  Wasliington,  DC,  Government  Printing  Office  (SD  Catalog  No.  13.10:348), 
1972,  p.  365-376. 
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Table  13.  X  ERROR  IN  DETERMINING  FLEXURE  STRESS 


a.  When  neglecting 
corner  radii 

b/d 

r/d 

1.0 

2.0 

4.0 

0 

0 

0 

0 

0.02 

0.1 

0.1 

0 

0.04 

0.4 

0.2 

0.1 

0.06 

0.9 

0.4 

0.2 

0.08 

1.5 

0.8 

0.4 

0.10 

2.4 

1.2 

0.6 

0.15 

5.1 

2.5 

1.3 

0.20 

8.6 

4.3 

2.2 

b.  When  neglecting 
45°  chamfer 

C/d 

1.0 

2.0 

4.0 

0 

0 

0 

0 

0.01 

0.1 

0.1 

0.1 

0.02 

0.2 

0.1 

0.1 

0.03 

0.5 

0.3 

0.1 

0.04 

0.9 

0.5 

0.2 

0.05 

1.4 

0.7 

0.4 

0.06 

2.0 

1.0 

0.5 

0.08 

3.4 

1.7 

0.9 

0.10 

5.2 

2.6 

1.3 

Note:  All  errors  are  negative. 


Load  Readout 

It  is  readily  apparent  that  an  error  in  the  break  load  P  is  identically  carried 

over  as  an  error  in  the  stress  o,  . 

b 

Specimen  Dimension  Measurement 

It  is  further  evident  that  an  error  in  measuring  the  specimen  dimension  can  also 
lead  to  an  error  in  stress.  It  is  recommended  that  the  cross  section  dimensions  b 
and  d  be  measured  at  the  point  of  failure  (to  preclude  specimen  taper  effects).  Con¬ 
sidering  the  true  specimen  dimension  to  be  in  error  by  e  ,  then  from  Equations  la  or 
lb:  m 


°_X _ bd* 

°b  ”  (b  +  eB)(d  7  eB)z 


for  three-  or  four-point  flexure. 


(24) 


If  em  is  small  relative  to  b  or  d: 

-  4  &)  *  ft)]  t2 

Equation  25  shows  that,  if  the  measurement  error  is  expressed  as  em/b  or  em/d,  the 
error  in  stress  is  magnified.  For  example,  if  d  *  b,  then  a  1%  error  in  specimen 
measurement  becomes  a  3%  error  in  stress. 
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RECOMMENDATIONS  FOR  FLEXURE  TESTING 


It  is  beyond  the  scope  of  this  report  to  analytically  determine  the  intersec¬ 
tion  of  errors  to  arrive  at  a  total  stress  error.  Nevertheless,  a  range  of  practi¬ 
cal  geometry  ratios  and  error  tolerances  can  be  specified  so  that  a  simple  additive 
stress  error  is  a  few  percent  at  most.  The  recommendations  are  summarized  in 
Table  14  and  are  discussed  below. 

Several  of  the  error  sources  are  negligible  for  most  common  test  configurations. 
These  include  initial  beam  curvature,  anticlastic  curvature,  beam  overhand  and  large 
deflection  sources.  The  error  due  to  nonhomogeneity  are  largely  unknown  at  th:*s  time. 

Many  of  the  errors  are  independent  of  the  test  configuration  but  should  not  be 
overlooked.  Micrometers  are  readily  available  that  are  accurate  to  within  0.0025  mm 
(0.0001"),  and  these  should  be  used  to  keep  specimen  dimension  measurement  errors  to 
a  few  tenths  of  a  percent.  Many  conventional  universal  testing  machines  can  easily 
read  break  load  to  within  0.5%.  Comer  chamfers  should  not  be  casually  applied  to 
specimens,  particularly  ones  with  small  cross  sections,  since  the  error  can  be  sig¬ 
nificant.  The  analysis  in  this  report  assumes  the  chamfers  were  identical.  If  they 
are  not,  or  if  only  two  chamfers  are  used,  a  further  error  can  result  due  to  a  shift 
in  the  position  of  the  specimen's  neutral  axis. 

Some  of  the  more  important  error  sources  do  depend  upon  the  fixture  configura¬ 
tion.  The  1/3-four-point  mode  has  somewhat  less  error  than  the  1/4-four-point  mode 
for  the  cases  of  wrong  span  and  contact  tangency  shift  sources.  A  greater  differ¬ 
ence  exists  for  the  eccentric  loading  source  of  error.  Special  care  should  be  taken 
to  minimize  wrong  spans  or  eccentric  loading  error  sources  in  four-point  flexure 
since  an  error  in  such  fixture  positioning  is  magnified  as  an  error  of  stress. 
Three-point  loading  is  much  less  sensitive  to  load  bearing  position  error  sources 
than  four-point  loading.  On  the  other  hand,  a  three-point  loaded  beam  is  adversely 
affected  by  the  presence  of  wedging  stresses  at  the  point  of  maximum  stress.  These 
wedging  stresses  decay  rapidly  with  distance  away  from  the  load  bearings  and  will 
have  considerably  less  influence  on  four-point  testing.  The  bearing  friction  error 
can  be  of  large  magnitude  for  either  three-  or  four-point  loading,  and  it  is  strong¬ 
ly  recommended  that  the  load  bearings  be  mounted  such  that  they  are  free  to  rotate. 
Twisting  error,  due  to  lack  of  parallelism  of  fixture  bearings  or  specimen  surfaces, 
is  harder  tc  predict,  because  the  error  is  dependent  upon  many  geometry  terms  as 
well  as  the  specimen  stiffness.  Parallelism  requirements  are  more  important  for 
four-point  loading  than  three-point.  For  most  geometries  and  materials,  parallelism 
limits  of  better  than  1°  in  the  specimen  and  also  the  fixtures  are  needed  to  keep 
the  error  within  1  percent. 

There  are  two  conflicting  requirements  regarding  contact  radius  at  the  loading 
and  support  points:  the  first  is  that  radii  must  be  great  enough  so  that  contact  or 
bearing  pressure  does  not  cause  local  failure  of  the  beam;  the  second  is  that  the 
contact  radii  be  small  enough  so  that  the  error  due  to  contact  point  tangency  shift 
is  not  great. 

Many  of  the  error  analyses  in  this  report  assumed  the  ratio  of  elastic  modulus 
to  bend  strength  (E/o)  was  1000.  Values  could,  in  fact,  range  from  100  to  2500.  In 
general,  the  larger  E/o,  the  larger  will  be  the  twisting  error  and  load  bearing  con¬ 
tact  stress,  but  the  lesser  will  be  the  contact  tangency  shift  and  large  deflection 
errors. 
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Table  14.  ERROR  SOURCES  AND  RECOMMENDATIONS 


Error  Source 


Recommendations 


Ef  t  E  ;  Nonhomogeneity  A  Anisotropy  (Table  1) 


Error  depends  upon  material  and 
on  a  fabrication  process. 


Initial  Beam  Curvature  (Table  2) 


Pc/d  >  100 


Antlclastlc  Curvature  (Table  3) 


b/d  <  T5 


Large  Deflection  (Table  4) 


Four-Point 

Three-Point 


a/d  <  12.5 
L/d  <  25 


Eccentric  Load  (Tables  5-7) 


1/3  Four- Point 

1/4  Four-Point 
Three- Point 


Non-Pi  voting  Head  e/L  <  0.001 
Pivoting  Head  e/’.  <  0.002 

Non-Pivoting  Head  e/L  <  0.001 
Pivoting  Head  e/L  <  0.002 

e/L  <  0.025 


Wrong  Span  (Table  8) 


1/3  Four- Point 
1/4  Four-Point 
Three- Point 


es/L  <  0.001 
es/L  <  0.001 
es/L  <  0.005 


Beam  Twisting  (Table  10) 


Minimise  es  and  ef  «  1° 


Bearing  Friction  (Eqs.  19,20) 


Four-Point 
Three-Pol nt 


Roller  bearings  which  are  free 
to  roll. 

Roller  bearings  for  outer  supports. 


Contact  Stress 


Four-Point 

Three-Point 


Pi/d  >  7.25  d/a 
Pi/d  >  29.0  d/L 


Wedging  Stress  (Table  11) 


Four-Pol nt 
Three- Point 


a/d  >  5.0 
L/d  >  20 


Beam  Overhang 


LT  >  L  A  2d 


Contact  Point 
Tangency  Shift  (Table  12) 


1/3  Four- Point 
1/4  Four-Point 
Three-Pol nt 


Pi/d  <  2.0 


Corner  Chamfer  (Table  13) 


Corner  Radius  (Table  13) 


Load  Readout 
Specimen  Dimension 


b/d  »  l.Oj  c/d  <.  0.03 

'b/d  »  2.0,  c/d  <  0.Q4 _ 

b/d  -  1.0,  r/d  i  0.04 
b/d  »  2.0,  r/d  <.0.06 
Measure  p  accurate  to  0.5X 
Measure  e^d  accurate  to  0.1X 


STRENGTH  AS  A  FUNCTION  OF  SPECIMEN  DIMENSIONS  AND  SAMPLE  SIZE 


General 

An  additional  issue  is  the  question  of  how  many  specimens  should  be  broken  in  a 
test  sequence.  Furthermore,  it  is  well  known  that  the  size  of  the  specimen  can 
influence  the  measured  strengths.  In  general,  the  larger  the  specimen,  the  weaker 
it  is  likely  to  be.  How  can  strength  results  generated  with  one  specimen  size  be 
compared  to  other  sizes?  These  two  questions  can  be  addressed  by  the  well-known 
Weibull  analysis.15 

Many  investigators  have  used  the  Weibull  approach  to  relate  strength  levels  of 
various  types  of  specimen  configurations  either  on  a  stressed  volume  or  surface  area 
basis.25  The  reader  is  cautioned  that  confirmation  of  such  an  analysis  or  lack 
thereof  may  well  depend  on  a  number  of  factors  including  the  test  material.  As 
examples  of  such  correlation  and  lack  of  it,  Weibull  statistical  correlation  was 
justified  by  Davies25  for  react ion- bonded  silicon  nitride  but  inappropriate  for 
Lewis'  work26  in  alumina  fabricated  by  several  processes. 

A  computer  program  for  statistical  evaluation  of  composite  materials,  applicable 
to  ceramic  materials,  is  available  in  Reference  27.  This  program  determines  the 
desirability  of  a  particular  probability  density  function  in  predicting  fracture 
strength  of  ranked  empirical  data.  The  candidate  functions  include  normal,  log 
normal,  and  Weibull.  Root  mean  square  error  results  can  be  tabulated  for  each 
functional  comparison.  The  effects  of  several  different  statistical  ranking  schemes 
can  be  readily  listed  in  the  computer  output. 

The  data  mean  and  standard  deviations  with  corresponding  levels  of  confidence 
can  be  included  in  the  printed  results.  The  Weibull  parameters,  obtained  from  the 
maximum  likelihood  method,  and  corresponding  confidence  intervals  can  be  obtained 
from  this  program. 

Since  a  Weibull-type  analysis  is  applicable  in  many  instances,  resulting  formu¬ 
las  for  the  simple  two-parameter  system25’28  to  determine  the  risk  of  rupture  for 
the  four-point  and  three-point  loading  systems,  are  presented  below,  for  the  sake  of 
completeness. 

It  is  worth  noting  that  Weibull  analyses  of  strength  data  require,  as  input, 
the  idealized  tensile  stress  acting  upon  a  specimen,  not  the  stress  at  the  point  of 
fracture.25  It  is  for  this  reason  that  strength  values  are  not  adjusted  in  four- 
point  loading  for  "out  of  inner  span  fractures"  (which  occasionally  occur)  or  for 
fracture  away  from  the  middle  bearing  in  three-point. 

Volume  Sensitive  Material 


The  Weibull  two-parameter  volume  distribution  function  for  the  probability  of 
failure  (F)  of  a  uniaxially  stressed  component  is:15*25 


25.  DAVIES,  D.  G.  S.  The  Statistical  Approach  to  Engineering  Design  in  Ceramics.  Proc.  Br.  Ceramics  Soc.,  no.  22, 1 973,  p.  429-452. 

26.  LEWIS,  D.,  III.  An  Experimental  Test  of  Weibull  Scaling  Theory.  J.  Amer.  Cer.  Soc.,  v.  59,  no.  1 1-12, 1976,  p.  507-510. 

27.  NEAL,  D.  M.,  and  SPRJGIDIGLIOZZI,  L  An  Efficient  Method  for  Determining  the  'A ’  and  ’B‘  Design  Allowables.  ARO  Report  83-2, 
Proc.  of  the  Twenty-Eighth  Conference  on  the  Design  of  Experiments  in  Army  Research,  Development  and  Testing,  1983. 

28.  DeSALVO,  G.  J.  Theory  and  Structural  Design  Applications  of  Weibull  Statistics.  Westinghouse  Astronudear  Laboratory,  WANL-TME- 
2688, 1970. 
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U  W  •] 


where  a  is  the  tensile  stress  acting  upon  an  element  dV  of  +he  component,  aQ  is  the 
characteristic  strength  (a  normalizing  parameter  which  has  units  of  stress  •  volume 
raised  to  1/M),  M  is  the  Weibull  modulus,  and  V  is  the  volume. of  the  component.  In 
general,  a  is  a  function  of  location  in  the  component.  Equation  (26)  is  often  re-  . 
written  for  flexure  specimens  in  terms  of  and  the  equivalent  volume  Vg  (the  vol¬ 
ume  of  a  tensile  specimen)  which,  when  subjected  to  the  same  stress  o^,  would  have 
the  same  probability  of  failure.25 


exp  -  - 


The  equivalent  volume  is  a  useful  quantity  since  it  permits  comparison  of  the 
mean  strengths  of  two  different  sized  components:* 


where  oi  and  V  .  are  the  strength  and  equivalent  volume  of  one  component,  and  o2  and 
Vg2  are  for  theother.** 

The  effective  volume  of  a  rectangular  beam  in  four -point  flexure  is: 


V  =  v  ■L.  -  . 
VE  2(M  +  1) 


,  2a  /  M 
1  “  L  M  +  1 


where  V  is  the  specimen  volume  inside  the  outer  span  (V  =  bdL) .  This  formulation 
includes  the  material  between  the  inner  and  outer  bearings.  For  the  case  of 
1/4-f our-point  bending: 


V  =  v  [  M+2 
E  4(M  +  l)2 


(30a) 


For  the  case  of  1/3-point  bending: 


vE  =  V 


r  m  +  3 

|_6(M  +  l)2 


(30b) 


•Or  the  itrengthi  »t  the  nine  probability  of  failure. 

••Some  asaumptions  are  involved  in  the  above  analysis;  the  leader  is  directed  to  Reference  25  for  details. 
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1 


and  for  a  bean  stressed  in  three-point  bending: 


»  vf - 1  , 

2(M  +  l)2 


(30c) 


M  can  be  determined  by  a  number  of  different  methods  (see  References  25,  28-30) . 

The  accuracy  by  which  M  can  be  determined  is  discussed  later  under  Weibull  Parameter 
Estimate  and  Sample  Size.  Equation  30  shows  that  a  larger  volume  of  material  is 
effectively  stressed  in  four-point  as  compared  to  three-point  loading.  It  is  for 
this  reason  that  four-point  loading  is  generally  preferred. 

Surface  Sensitive  Material 

The  Weibull  two-parameter  surface  distribution  function  for  the  probability  of 
failure  of  a  uniaxial ly  stressed  component  is:15*25 


i  -  «p  -  r  ~ 


(31) 


where  the  characteristic  strength  has  units  of  (stress  •  area  raised  to  1/M),  M  is 
the  Weibull  modulus,  and  integration  is  performed  over  the  specimen  surface,  S.  If 
surface  flaws  predominate,  then  the  effective  surface  Sg  can  be  used  to  compare 
mean  strengths*  between  two  components: 


2.  .  I  f|2\  1/M 

02  '  \SE1/ 

Tiu  Tfective  surface  area  for  a  four-point  loaded  beam  is: 

2 

SE  “  (m4t)  |bA(M  +  1)2  +  C2ab  +  dAl(M  +  D  +  2adj  » 

and  for  three-point  it  is: 

„  Lb  Ld 


(32) 


(33) 


E  M+1  (M  +  l)2 

Once  again,  a  greater  surface  is  exposed  to  high  stress  in  four-point  loading  as 
compared  to  three-point,  which  is  why  four-point  is  preferred. 


(34) 


*Or  the  strengths  at  the  tame  probability  of  failure. 

29.  McLEAN,  A,  F.,  and  FISHER,  R  A.  Brittle  Materials  Design,  High  Temperature  Gas  Turbine.  Ford  Motor  Company,  Contract 
DAAG46-7 1-00162,  Interim  Report,  AMMRC  CTR  77-20,  August  1977. 

30.  PALUSZNY,  A.,  and  WU,  W.  Probabilistic  Aspects  of  Designing  with  Ceramics.  Presented  at  the  22nd  Annual  Gas  Turbine  Conference 
of  A.S.M.R,  Philadelphia,  Pa.,  March  27-31, 1977. 
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Mel bull  Parameter  Estimate  and  Sample  Size 


Often,  the  objective  of  flexural  strength  testing  is  not  merely  to  estimate  a 
mean  strength,  but  to  estimate  Weibull  distribution  parameters  such  as  the  Weibull 
modulus  M,  or  the  characteristic  strength,  o  .  The  following  section  discusses 
requirements  for  numbers  of  specimens  in  order  to  obtain  reasonable  estimates  for 
these  parameters. 

Flexure  tests  on  hot-pressed  silicon  nitride  material  reported  by  McLean  and 
Baker31  show  the  effect  of  Weibull  slope  M  for  specific  component  reliability.  The 
strength  requirement  for  a  specific  component  reliability  was  decreased  16%  by  a 
reported  20%  increase  in  M  from  a  nominal  value  of  10,  and  was  increased  27%  by  a 
20%  decrease  in  the  slope. 

Different  techniques  will  produce  somewhat  different  results,  according  to 
McLean  and  Fisher,29  when  estimating  the  Weibull  parameters.  Two  statistical  methods 
had  been  used  during  preliminary  analysis  of  hot-pressed  silicon  nitride  material 
strength  data,  and  results  indicated  that  the  estimates  of  the  characteristic  value 
ap  (or  scale  parameter)  were  very  close  while  the  Weibull  slope  estimates  vary  and 
thus  would  yield  considerable  differences  in  the  component  strength  requirement . 

The  following  is  quoted  directly  from  Reference  29  (except  to  change  reference 
and  figure  numbers  approprirte  for  this  report)  because  it  succinctly  addresses  the 
answer  to  the  question  of  proper  sample  size:  'The  exact  confidence  intervals  for 
the  parameters  are  based  on  the  distributions  obtained  by  Monte  Carlo  methods  presented 
in  Thoman  et  al.32  It  is  not  unexpected  that  the  uncertainty  in  the  estimation  of  a 
parameter  will  increase  as  the  sample  size  decreases.  This  uncertainty,  however,  has 
rarely  been  quantified.  The  width  of  the  confidence  intervals  for  the  parameters  is 
a  measure  of  the  uncertainty  and  aids  in  the  selection  of  the  sample  size  of  a  test. 
Figures  6  and  7  are  drawn  from  Reference  32  and  show  the  90%  confidence  bounds  for 
the  Weibull  slope  and  the  characteristic  value."  (Figure  7  differs  from  that  given 
in  Reference  29  in  that  two  additional  M  values  were  computed  and  shown.)  "The  bounds 
for  the  Weibull  slope  are  a  function  of  sample  size  only,  while  for  the  characteristic 
value  they  are  a  function  of  both  the  sample  size  and  the  Weibull  slope.  As  can  be 
seen  from  the  graphs,  the  error  or  uncertainty  in  estimates  from  small  sample  sizes 
is  very  large.  Important  judgements  and  significant  analysis  should  not  be  based  on 
small  samples.  Sample  sizes  of  at  least  30  should  be  used  for  all  but  the  most  pre¬ 
liminary  investigations.  An  uncertainty  of  ±10%  in  Weibull  slope  requires  more  than 
120  samples.  This  uncertainty  is  not  peculiar  to  just  ceramics,  but  is  intrinsic  to 
the  statistical  analysis  of  data,  whether  that  data  be  material  strength  or  the  life 
of  some  electronic  component.  The  choice  of  sample  size  depends  on  many  factors  in¬ 
cluding  the  cost  and  timing  of  testing  and  the  degree  of  conservation  which  is  accept¬ 
able,  but  erroneous  judgments  may  be  made  and  unacceptable  designs  pursued  if  the 
sample  sizes  are  too  small." 


31.  McLEAN,  A.  F.,  and  BAKER,  R.  R.  Brittle  Material!  Design,  High  Temperature  Gas  Turbine.  Ford  Motor  Company,  Contract 
DAAG46-7 1-00 1 62,  Interim  Report,  AMMRC  CTR  76-31,  October  1976. 

32.  THOMAN,  D.  R.,  BAIN,  L.  J.,  and  ANTLE,  C  E  Inferences  on  the  Parameters  of  Weibull  Distribution.  Technometrica,  v.  11,  p.  445-460. 
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Characteristic  Value  Error  Percent  of  Expected  Value 

•  *-*  02 


Figure 7*  aw«rei«lc»«lu>««^*^to-TOo°n,,^1*n*' 
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LOADING  SPEED 


It  is  well  known  that  speed  of  loeding  will  often  influence  the  failure  stress 
of  structural  ceramic  beams.  The  source  of  the  sensitivity  is  stress  corrosion 
phenomena,  particularly  in  the  presence  of  water  or  water  vapor.  In  general*  the 
slower  the  speed  of  loading,  the  greater  the  opportunity  for  stress  corrosion 
phenomena  to  weaken  the  specimen.  Thus,  fast  loading  speeds  are  usually  used  in 
strength  tests. 

Most  universal  testing  machines  used  for  flexure  testing  are  constant  displace¬ 
ment  rate  machines,  so  it  is  convenient  to  specify  strain  rates.  The  strain  rate 
for  a  linearly  elastic  material  is  defined  as: 

e  •  (ob/E)/t  ,  (35) 

where  t  is  the  time  of  the  applied  load;  but  since  the  speed  of  loading  is  s  *>  y/t, 
then 

i-  (£)(*).  CM) 

where  y  is  the  deflection  of  the  beam  and  s  the  constant  speed  of  the  testing 
machine.  This  assumes  that  all  of  the  machine's  motion  is  transmitted  to  the  spec¬ 
imen  (i.e.,  the  machine  is  perfectly  'hard').  The  deflection  of  the  inner  load 
bearings  of  a  four-point  loaded  beam  is: 


y  *  (3L  -  4a)  . 

Substitution  of  y  into  Equation  36  and  recognizing  that  ob  ■  Pad/21  gives: 

•  3Sd 

e  "  a(3L  -  4a)  * 

For  a  1/4- four-point  beam: 

•  m  6dS 
€  "  L2 

For  a  1/3-point  beam: 

•  27dS 

e  »  -  . 

5L2 


(37) 


(38) 


(39) 


(40) 


Using  the  same  approach  as  above,  we  obtain  the  strain  rate  for  three-point  beams: 


6dS 
e  »  - 

L2 


(41) 
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An  alternative  approach  to  designating  the  loading  speed  is  the  stressing  rate: 

o  -  cB  (42) 

which  is  valid  for  the  oase  of  a  linearly  elastic  material.  Equations  38  to  41  can 
be  substituted  into  Equation  42  if  stressing  rates  are  specified* 

Finally,  if  the  tiae  per  test  is  the  limiting  concern,  then  the  following  is 
applicable: 

t-'T*  (d3) 

and  Equations  38-42  can  be  used  along  with  the  projected  bend  strength  to  solve  for 
t. 


CONCLUSIONS 

A  variety  of  sources  can  lead  to  errors  in  determining  the  flexure  strength 
when  using  simple  beam  theory  equations.  These  sources  include  assumptions  involv¬ 
ing  simple  beam  theory  and  external  influences  pertaining  to  the  load  application. 
Providing  that  the  beam  is  homogeneous  and  isotropic,  and  deflections  are  relatively 
small,  then  the  major  sources  of  error  are  from  external  sources.  In  particular, 
the  most  serious  errors  arise  from  load  bearing  friction,  beam  twisting,  and  load 
bearing  mislocation.  Other  errors,  such  as  contact  point  tangency  shift,  wedging 
stresses,  neglecting  comer  chamfers,  and  load  readout  errors  cannot  be  neglected 
either.  Table  14  lists  all  of  the  potential  error  sources  identified  in  this  report 
and -makes  speoific  recommendations  for  specimen  and  fixture  geometries  and  toler¬ 
ances.  The  bases  for  the  recommendations  are  that  they  be  practical,  that  they 
limit  the  individual  errors  to  approximately  one  half  percent,  and  that  the  sum  of 
the  errors  be  less  than  a  few  percent. 

Requirements  for  a  minimum  number  of  recommended  specimens  (30)  are  presented 
in  the  context  of  the  Weibull  two-parameter  analysis.  This  analysis  is  one  of  the 
simplest  possible,  and  the  reader  is  cautioned  that  numerous  assumptions  are  entailed 
in  its  use.  Even  if  a  more  complex  function  appears  to  have  better  applicability 
than  a  Weibull  analysis,  the  requirement  for  30  or  more  specimens  should  likely  re¬ 
main  valid. 

For  convenience,  a  brief  discussion  of  converting  strength  of  one  size  speci¬ 
men  to  another  is  included.  Again,  since  this  analysis  is  based  upon  a  Weibull  two- 
parameter  approach,  the  reader  is  cautioned  that  numerous  assumptions  apply  and  that 
more  sophisticated  analyses  may  have  to  be  used. 

A  section  on  loading  speed  is  also  included  for  convenience  to  permit  quick 
assessment  of  optimum  universal  testing  machine  speeds. 
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TABULATIONS  OF  ERRORS  IN  CALCULATING  FLEXURE  STRESS 


_  Unless  otherwise  stated,  the  percent  error  is  determined  throughout  the  text  as 
*  *  t(0.  -  oj/o  1100;  where  o.  ■  6M/bd2  and  o  is  more  nearly  the  true  bending 
stress.  Thus  a  negative  error°indi cates  the  simple  beam  formulas  la  and  lb  under¬ 
estimate  the  true  stress;  a  positive  error  is  an  overestimate. 
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APPENDIX  A.  ANTICIASTIC  CURVATURE 

When  a  ban  is  bant  by  a  moment,  it  produces  a  curvature  p  along  its  longitu¬ 
dinal  axis;  there  is  also  curvatura  prasant  in  tha  transvarsa  or  lataral  direc¬ 
tion.  This  nonant  is  definad  by  orthodox  theory  as 

%  -  (BI/p)6  (A-l) 

where  Ms  a  paranater  representing  tha  affect  of  restraint  of  anticlastic  curvatura 
after  Ashwell.9  Since 


°x  “  N^c/I  *  (BI/pI)y& 

than 


ox  *  (By/P)3  (A-2) 

where  y  is  the  distance  from  the  neutral  axis. 


The  8  for  siuple  bean  theory  will  equal  1.0  and  if  the  bean  width  to  depth  is 
great »  i.e.,  b/d*+** ,  the  bean  can  be  considered  as  a  plate  so  that  B*l/(l-v*).  It 
is  worthwhile  to  know  the  intermediate  values  of  g  such  that  the  effect  of  restraint 
of  anticlastic  curvature  on  the  error  can  be  ascertained  when  assuning  simple  bean 
theory  (8  »  1.0)  is  valid. 


Ashwell9  has  determined  g  as  a  function  of  Poisson's  ratio,  bean  width, 
depth,  and  neutral  axis  curvature  by  accounting  for  anticlastic  curvature  and 
treating  the  structure  as  a  beam  on  an  elastic  foundation.  The  function  6  and 
related  terns  are  repeated  here  in  the  following: 


8 


where 


y 


2  yb 


f(Yb) 


2/3  v 
yb/l-v^ 


F(yb) 


f(yb)  ■  2(B2+C2) [sinh(yb)  ♦  sin(yb)] 

♦  (B2-C2*2BC)  cosh (yb)  sin(yb) 

♦  (B2-C2-2BC)  sinh(yb)  cos (yb) 

♦  2(B2-C2) (yb) , 

F(yb)  ■  (B+C)  sinh (yb/2)  cos(yb/2) 

-  (B-C)  cosh (yb/2)  sin (yb/2). 


(A-3) 


•To  be  non  met,  b2/pd  ♦  - . 
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By  programming  (A-3)  and  prescribing  b/d,  but  first  allowing  6*1.0,  then 
iterating  in  the  computer  through  (A-4),  the  relationship  between  b/d  and  B  was 
obtained.  Once  this  relationship  is  known,  the  percent  error,  defined  as 
[ (1-B)/B] 100,  as  a  function  of  b/d  is  realized.  These  errors  are  given  in 
Table  3  as  a  function  of  b/d  with  E/a j,  *  1000  and  v  *  0.25. 


APPENDIX  B.  LOAD  MISL0CATI0N  ERROR,  LOADING  HEAD  RIGIDLY  ATTACHED 
Four-Point  Loading 


Consider  the  usual  flexure  testing  setup  where  the  loading  head  is  rigidly 
attached  to  a  testing  machine,  schematically  shown  in  Figure  lb  and  idealized  in 
Figure  lc.  The  upper  loading  head,  where  the  inner  span  t  is  fixed,  can  only 
translate  in  the  vertical  direction,  and  the  lower  support  fixture,  where  the 
outer  span  L  is  fixed,  can  be  located  with  reference  to  the  loading  head.  The 
slope  and  deflection  equations  between  points  AB,  BC,  and  CD  are  as  follows: 


M  (dyAB/dx)  -  (Pix2/2)  ♦  c^ 
El  yAB  *  (PiX3/6)  +  cix  +  c2< 


(B-la) 


El  (dyi  lx)  »  (P1-P2) (x2/2)  +  P2aix  +  c3 
El  YBC  “  (pl-p2)(x3/6)  +  (P2aiX2/2)  +  c3x  +  c4 


(B-lb) 


and 


El  (dyco/dx)  -  (pi-P2-p3)(x2/2)  ♦  P2aix  +  P3  (ai *t)  x  +  c5 

BI  yco  -  (pi-p2-p3) (x3/6)  +  (p2»ix2/2)  +  P3  (ai+£) (x2/2)  +  c5x  ♦  c6 


(B-lc) 
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where  Pi,  P2,  P3,  and  P4  are  the  loads,  x  and  y  are  defined  as  shown  in  Figure  1, 
E  is  the  Young's  modulus  of  the  material,  and  1  is  the  moment  of  inertia  of  the 
cross  section  of  the  beau. 

Through  the  use  of  the  various  boundary  conditions  the  constants  are  deter¬ 
mined  to  be: 

d  «  -1/6  {Pi (L2-aj)  +  (P2/L)  [(ai-L)3  +  (f+a2)a22]} 
c2  *  0 

c3  -  -1/6  {Pi(L2-a22)  +  (P2/L)  [3ai2L  -  (a 2+l)3  +  (£+a2)a22]} 

C4  =  P2ai3/6 

c5  =  -1/6  {PXL2  +  (P2/L)  [ai3+3aiL2-L3]  +  (P3/L)  [(a^*)3  +  3(ai+£)  L2-L3]> 
c6  »  P3[(ai+^)3/6]  +  (P2ai3/6). 

In  order  to  determine  the  distribution  of  loading  between  the  vertical  loads 
Pi,  p2»  p3»  and  P4,  the  final  condition  of  equal  deflection  must  be  enforced  at 
locations  B  and  C  (see  Figure  lb) ,  which  is 


(yBC)x«ai  "  ^■yBC)XtIai+^< 

Enforcing  the  above  condition  in  the  second  part  of  Equation  B-lb  results  in 


pl/P2  -  -tf/L)2 


If L  -  (1-ai/L) (2-£/L-2ai/L) 

(ai/L)3  -  (ai/L+^/L)3  ♦  (1  -  [l-(A/L)-ai/Li2)(£/L)  . 


(B-2) 


Utilizing  force  and  moment  equilibrium,  a  further  relationship  between  all  four 
forces  is  obtained  and  given  in  the  following: 

CPi/p2)(ai/L)  ♦  (Pi/P2-l)f/L  3  ^  ^  Px/Pz  +  ai/L  -  1 

P4/P3  ■  - ,n  /n  - - 7— - y. -  and  P3/P2  » -  (B-3) 

(P1/P2)  "  1  +  ai/L  1  -  ai/L  -  t/L  . 

The  ratio  of  the  stress  at  x  (ax)  to  the  bending  stress  (ob)  from  Reference  16 
or  Equation  9  in  the  text,  where  it  is  assumed  that  ai  /  a2,  is: 


where  a  is  the  value  ut  ai  With  perfect  load  location  and  xi  is  defined 
as  shown  in  Figure  lc. 

By  manipulation  of  (B-2)  and  (B-3),  the  factor  Pi/(P2+  Ps)  in  (B-4)  can  be  put 
into  terms  of  ai/u  and  l/L.  This  was  accomplished  and  the  results  are shown  in  the 
following  equation: 
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V°b 


[  (ax/L)/(a/L)]  (1  -  ax/L  -  A/L) [(A/L) (2  -  a1/L)-2(l  -  ax/L)2] 
3(ai/L)(l  -  A/L  -  ax/L)  -  (1  -  A/L)2 


(B-5) 


By  defining  the  eccentricity  of  loading  as  e/L  *  ax/L  -  a/L,  Equation  B-5  becomes: 

[(e/L  +  a/L)/ (a/L)] [1- (e/L  +  a/L)-  A/L] { (A/L) [ 2- (e/L  +  a/L)]- 2[l-(e/L  +  a/L)]2) 

ax/°b  * - — - - —  _ : _ 

3 (e/L  +  a/L)  [1  -  A/L  -  (e/L  +  a/L)]-  (l-  A/L)2 


(B-6) 


Calculations  of  ox/ob  were  obtained  for  A/L  =  1/3  as  well  as  A/L  =  1/2,  by 
allowing  e/L  to  take  on  negative  values  only  in  Equation  B-6.  Only  negative  values 
were  considered  because  beam  failure  will  occur  due  to  a  realistically  larger  moment 
than  idealized  when  ignoring  eccentricity.  These  error  calculations,  although  de¬ 
termined  by  allowing  e/L  <  0  in  Equation  B-6,  are  indicated  as  ±e/L  in  Table  5 
for  A/L  =  1/3  and  Table  6  for  A/L  *  1/2.  This  simply  indicates  that  the  location  of  the 
maximum  moment  or  stress  is  at  xx  =  ax  when  e/L  <0  and  at  xx»  ax  +  A  when  e/L  >  0. 


The  reader  is  cautioned  that  for  each  value  of  a/L  there  exists  a  set  of  limits  on 
(B-2 ),  (B-3),  and  (B-5).  That  is,  ax  can  be  such  that  either  P2  or  P3  can  equal  zero, 
because  O^bC^x^  /  O^BC^xaai+A  and  t^e  system  changes  from  four-point  to  an  eccentric 
three-point  loading.  The  limiting  values  can  be  determined  by  allowing  P2  =  0  in  Eaua- 
tion  B-2. 


APPENDIX  C.  BEAM  TWISTING 

If  line  loads  are  nonuniform  or  nonparallel  between  pairs  of  load  contacts, 
or  if  the  cross  section  of  the  specimen  is  skewed  along  its  length,  as  shown  in 
Figure  3,  a  net  torque  will  result.  The  addition  of  torque  gives  rise  to  a  max¬ 
imum  principal  stress  due  to  bending  and  torsional  stresses. l3>*6  Failure  assumed 
to  be  caused  only  by  bending  stress  will  yield  an  error.  Two  cases  are  consi¬ 
dered:  Case  I  -  Failure  occurs  prior  to  specimen  realignment  in  the  bend  fixture 
(bottoming),  and  Case  II  -  Failure  occurs  at  or  after  bottoming. 

Case  I 

Recalling  that  the  bending  stress  for  a  loaded  beam  is 

ax  =  ab  =  6Mb/bd2  (C-l) 

where  is  the  measured  bending  moment  at  failure;  for  a  four-point  loaded  beam 
Mjj  *  Pa,  and  for  a  three-point  loaded  beam  =  PL/4. 

The  maximum  shear  stress  due  to  torsion  of  a  rectangular  beam  is3 

Txz  "  Vk2bd2  Cc-2) 

where  Tb  is  the  torque  and  equal  to  Pb  for  four-;,oint  bending  and  Pb/2  for  three- 
point  bending,  and  kg  is  a  numerical  factor  obtained  from  Reference  3  and  is  given 
in  Table  9.  This  peak  shear  stress  occurs  at  the  specimen  surface  at  the  midpoint 
of  the  long  edge  (dimension  b). 
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Prior  to  bottoming,  th«  normal  stress  is 


(g^Hl-cosie) 


♦  (0^/2}  (l+cos28)  ♦  txl#in26 , 


(C-3) 


on  a  plane  whose  normal  is  in  the  xz  plane  and  is  inclined  at  ah  angle  6  to  the  x 
axis.  Since  we  shall  assume  a  plane  strain  condition  i.e. ,  e  •  0,  then 

m 

e  *  0  ■  (1/E) (o  -vb  )  or  ow  ■  vo  ■  vc 
since  cr^  «  0  at  the  free  surface.  From  Equation  C-3: 

°n  *  (ob/2) [(1+v)  +  (l-v)cos20]  +  Txzsin20.  (C-4) 

Now  a  is  maximum  when 
n 

tan26*  «  2tXz^°x“°z^1 
but  since 


then 


tan20*  -  2txz/ (!-v)Ob  (C-5) 

where  6*  is  the  angle  of  a  plane  inclined  to  the  axis  at  which  principal  stress  is 
a  maximum.  Substitution  of  (C-l)  and  (C-2)  into  (C-5)  for  the  condition  prior  to 
bottoming  gives: 

tan20*  -  (Tb/Mb)/[3(l-v)k2L  (C-6) 

sin20*  -  (V^/lCy^)2  ♦  (3k2)2(l-v)2]1/2,  (C-7) 

and 

cos20*  -  dMl-v^/tCy^)2  +  (3k2)2(l-v)2]1/2.  (C-8) 

From  Equations  C-l  and  C-2 

Txz  "  (V2)[(Tb/Mb)/3k2l»  (C_9) 

and  by  substitution  of  the  above  relationships  into  (C-4)  with  some  algebraic  manip¬ 
ulation,  we  obtain: 

\ax  "  (°b/2)  iU+V)  +  (1/Sk2)tCW2  +  9C1-v)21c2231/2}  (CIO) 

prior  to  bottoming.  The  shear  stress  due  to  torsion  can  be  related  to  the  twist 
angle3  of  the  beam  through  the  following  relationship: 

*xz  -  CklAzJGdK^/L,.)  ♦  dp/*’)!  CC-U) 
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where  $s  is  the  twist »angle  along  the  length  of  the  specimen  (see  Figure  3b),  fp  is 
the  twist  angle  between  a  pair  of  load  and  contact  points  relative  co  $s  (see  Figure 
3c),  V  is  equal  to  either  "a"  for  four-point  beam  systems  or  L/2  for  three-point 
beam  systems,  ki  is  another  numerical  factor3  given  in  Table  9,  and  G  »  E/2(l*v), 
the  shear  modulus  of  the  material. 

Equation  (C-2)  can  be  equated  to  (C-ll)  and  thus  we  obtain: 

Tb  -  [k1E/2(l+v)]bd2[(d/LT)^s  +  (d/Jt»)*p]  (C-12) 

e 

where  T.  is  the  torque  when  bottoming  occurs.  Thus,  in  order  for  (C-10)  to  be 
De 

applicable,  T.  must  be  less  than  T.  ,  and  since  T,  «  (b/i')M.  and  from  (C-l)  and 
(C-12):  b  be  bo 

Tb  /Tfe  »  [3kiE/ab(l+v)(b/f)]bd2[(d/LT)^s  +  (d/t'Hp]  >  1.0,  or 


n«  (JlVb)[3k1(E/ab)/(l+v)3[(d/L..)*s  +  (d/f)*p]  >  1.0 
where  ob  is  the  bending  stress  at  failure  according  to  (C-l) . 

Note  that  T, /M,  *  b /£'  and  for  four- point  loading  Jt,'  =  a;  for  three-point 
loading  V  =  L/2,  tnus  (C-10)  becomes: 

»  (ob/2)  j(l+v)  +  (l/3k2)[(b/4*)2  +  9(l-v)2k22]1/2J 
with  n  >  1.0. 


(C-13) 


(C-14) 


Case  II 


If,  however,  n  <  1.0  then  bottoming  occurs  prior  to  or  at  failure  and  the 
following  analysis  is  applicable. 

Equations  (C-4)  and  (C-5)  are  still  appropriate  but  the  shear  stress  is 

Txz  *  Tb  /k2bd2*  (C-1S) 

e 


Substitution  of  t  from  the  above  and  ob  from  (C-l)  into  (C-5)  gives 

tan26*  -  (Tfe  /Mb)/3(l-v)k2, 

© 

but  from  (C-13)  T,  =  nT.  and  thus: 
be  b 

tan26*  -  (nTb/Mb)/3(l-v)k2.  (C-16) 


37 


As  in  Case  I,  using  a  like  procedure  we  determine  a  to  be: 

....  "max  .... 

°n  ~  (ob/2)  {  (1+v)  +  (l/SkajEn2^/^)2  ♦  9(l-v)2k22l1/2  j  (C-17) 


with 

n  “  (t*)  [3kl(E/ob)/C1+v)1[(d/LT^s  +  (d/A,)*P]  1  1-°‘  (C-18) 

Equation  (C-17)  is  applicable  to  both  systems  since  (Th/M^)  “  b/t'  *  for  the 
four-point  beam  bending  system,  and  (T.  /M. )  *  2b/L  for  the  three-point  system.  Thus 
(C-17)  becomes  d  d 

°n  B  (°b/2)  !(1+v)  +  U/3k2)[(nb/Jl*)2  +  9(l-v)2k22]1/2J  .  (C-19) 

max  *  ' 

Notice  when  n  =  1,  (C-19)  reduces  to  (C-14). 

Finally,  the  percent  error  is  defined  as: 

7  *  [(o.-o  )/o  ]100.  (C-20) 

max  max 

Errors  were  calculated  in  accordance  with  (C-20),  with  v  *  0.25  for  Case  I  (n  *  0.20, 
0.40,  0.60,  and  0.80)  and  Case  II  (n  ®  1.0). 


If,  instead  of  plane  strain  (ez  3  0),  it  had  been  assumed  that  a  plane  stress 
(oz  s  0)  condition  applied,  then  the  maximum  principal  stresses  are  given  by  the 
same  Equations  C-10,  C-14,  C-17,  and  C-19,  but  with  Poisson's  ration  v  »  0.  Equa¬ 
tions  C-12,  C-13,  and  C-18  are  unchanged  however,  since  G  ■  E/2(l+v).  The  plane 
strain  analyses  gives  a  higher  error  estimate,  but  the  plane  stress  condition  is 
closer  to  the  actual  case  since  lateral  constraint  is  negligible. 


APPENDIX  D.  WEDGING  STRESSES 

We  allow  the  stress  in  the  x  direction  in  Figures  1  and  2  to  be 

°x  a  °b  +  (2P/bd)BT,  (D-l) 

where  ob  is  the  bending  stress,  i.e.,  ab  »  (6Mx/bd2),  and  2P/bd  is  the  local 
stress,  i.e.,  the  so-called  wedging  stress,  and  8x  is  a  numerical  factor  dependent 
on  the  normalized  distance  x'/(d/2)  on  either  side  of  the  applied  load  point. 2 

For  convenience  the  value  of  B-r  at  the  tensile  side  of  the  beam  as  a  function 
of  x'/d  is  given  in  Table  D-l. 

The  percent  error  is  defined  as: 

[(<?b-ox)/ox]100,  (D-2) 
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and  substitution  of  (D-1)  into  the  above  equation  gives 

r-  ^-frr/Kob^p/bdj+^l^ioo.  (d-3) 

Since  Ob  *  6Mx/bd2,  then  (D-3)  becomes: 


7  .  ^-eT/(3Mx/Pd+8x)|  100.  (D-4) 

For  a  four-point  loaded  beam  the  bending  moment  is  constant,  i.e.,  Mx  ■  Pa, 
and  thus  equation  (D-4)  becomes: 

e  »  {-0T/  [  (3a/d)  +0T]}lOO.  (D-S) 

From  Table  D-1  and  Equation  D-5  above,  it  is  seen  that  the  error  is  dependent 
on  0x  or  the  fracture  location,  which  is  the  normalized  distance  x'/d.  These 
errors  have  been  computed  for  the  four-point  loaded  beam  and  presetted  in  Table  11a. 

For  the  three-point  loaded  beam,  (D-4)  is  still  applicable,  but  recalling  that 
Mx  »  P/2[(L/2)-x']  and  substituting  Mx  into  (D-4)  gives: 


=  ^-px/(l(L/d)-|(x’/d)+8T)j-100  for  x*  <  L/2. 


(D-6) 


Again,  as  can  be  seen  by  (D-6),  the  percent  error  is  dependent  upon  0T  and 
the  normalized  fracture  location.  These  percent  errors  are  given  in  Table  lib. 


Table  D-1.  WEDGING  PARAMETER  pT 


±x'/d  eT 


-0.1332 

+0.0137 

+0.0868 

+0.0640 

+0.0421 

+0.0220 

+0.0095 

+0.00075 


APPENDIX  E.  CONTACT  POINT  TANGENCY  SHIFT  , 

Consider  the  four-point  loaded  beam  shown  in  Figure  4.  The  original  span  "a" 
is  seen  to  decrease  by  the  amount  (hi+h2)  due  to  rolling  or  slipping  of  the  beam 
on  its  support  and  load  points.  The  beam  fulfills  the  condition:* 


0 

0.125 

0.250 

0.375 

0.500 

0.750 

1.000 

1.500 


*It  to  readied  that  a  chanfe  in  the  neuttant  load  vector*  mint  live  itoe  to  horizontal  reaction*;  hut  it  to  amumed  that  till*  affect  to  amalL 
and  to  thua  ignored. 
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dyVdx2  ■  Mx/E  I . 

The  moments  ere  defined  es: 

Mx  -  P(x-hi),  0  <  x  S  (a-h2) 

M*  ■  P[a-(hx+h2)],  (a-h2)  <  x  <  L-(a-h2). 

The, slope  equations  are 

El (dy/dx)  -  Pt(x2/2)-hxx]  +  Cj,  0  <  x  S  (a-hx),  and 

EI(dy/dx)  -  P[a- (hi+h2)]x  +  C2,  (a-h2)  <  x  <  L-(a-h2). 

Now  when  x  =  a-h2, 

cl  ♦  P  ^[(a“h2)2/2]  -  (a-h2)hx^  *  P[a-(hi+h2) J(a-h2)  +  C2,  or 
C2-Cx  -  “(P/2) (a-h2)2. 


(E-l) 


(E-2) 

(E-3) 


(E-4) 


Note  also  that  when  x  *  L/2  and  dy/dx  *  0  in  (E-3), 

C2  0  -P[a-(hx+h2)]L/2. 

After  substitution  of  C2  into  (E-4)  we  obtain: 

Cx  -  (P/2)  |(a-h2)2  -  L[a-(hi+h2)]J.. 

Substitution  of  th-se  constants  into  the  appropriate  slope  equations  gives: 

El  (dy/dx)  -  P[(x2/2)-hix]  +  (P/2)  {(a-h2)2  -  L[a-(hj+h2)] }  (E-7) 

with  0  <  x  <  (a-h2),  and 


(E-5) 

(E-6) 


El (dy/dx)  -  P[a-(hi+h2)][x-(L/2)], 


(E-8) 


with  (a-h2)  <  x  <  L-(a-h2). 

Now  when  x  *  ,-j  frr;  _  e  geometry  dy/dx  *  -hi//Pi-h2  •  -hi/Pi  and  when 
x  ■  a-h*,  dy/dx  -  -ha/  <  2-h|  <•  -ha/P2.  The  above  relationships  are  used  with  (E-7) 
and  (E-8)  and  we  obtain: 


(Pl/d) 


(hx/a)  (E/ob) _ 

(hx/a)2  -  (  *  a)2  ♦  (L/a)[l-(hx/a  +  h2/a)] 


(E-9) 
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•ad  .  ■  •. 

l-(h2/a)  •  (l/2){[(L/2s)  ♦  (hi/a)  ♦  A2) 

-  /[(L/2*)  ♦  (hi/a)  ♦  A2}2  -  4[Ch1/*)(L/2a)  ♦  A2]}  (B-40) 

Where  A2  ■  (E/2ob)/(p2/d). 

Note  that  for  the  thi se-point  loading  case  h2  ■  0,  L/a  »  2,  and  P  P/2  and 
(E-9)  reduces  to: 

;  ■  ■  i  ■  *  ' 

Pl/d  -  [2(h1/L)CB/ob)]/((2h1/L)-l]2,  (B-ll) 

and  the  region  of  validity  of  (E-10)  vanishes. 

The  percent  error  is  defined  as: 

e  -  [Cob-ax)/ox]100  -  [(^-M*) /Mr]  100,  or 

e  -{[(hj/a)  ♦  (h2/a)]/[l  -  (hi/a)  -  (h2/a)]}l00  (B-12) 

for  four-point  loading  and 

7  »{(2hi/L)/[l  -  (2h1/L)]}lOO  (B-13) 

for  three -point  loading. 

Calculations  were  performed  by  the  following  procedure.  It  was  assumed  that 
E/ob  ■  1  x  105,  and  thus  A*  ■  (1  x  10s)/(p2/d);  then  for  the  1/3-  and  1/4-four- 
point  loading  case  a/L  ■  1/3  and  1/4,  numerical  values  were  assigned  to  p2/d  and 
hi/a  and  then  h2/a  was  determined  from  (B-10).  Those  numerical  values  at  hj,/a 
and  corresponding  h2/a  were  substituted  into  (B-9)  to  determine  pi/d.  This  same 
procedure  was  used  for  the  three-point  loading  case  with  h2  «  0  and  L/a  «  2.  Once 
the  parameters  hi/a,  h2/a,  p2/d  and  Pi/d  are  known,  percent  errors  according  to 
(E-12)  for  the  four-point  loading  case  and  (E-13)  for  the  three-point  loading  case 
can  be  determined.  Such  errors  are  given  in  Tables  12a-c. 


APPENDIX  F.  ERROR  DUE  TO  NEGLECTING  CHANGE  IN  MOMENT  OF  INERTIA 
CAUSED  BY  CORNER  RADII  OR  CHAMFERS 

Comer  Radius 

Consider  Figure  5a,  which  shows  the  cross  section  of  a  rectangular  bean  with 
comer  radii  r.  The  true  moment  of  inertia  (I  „)_  about  the  centroidal  or  neutral 
axis  x-x  is: 

(I  -  b(d-2r)s/12  ♦  (b-2r)rs/6  ♦  (l/2)(b-2r)(d-r)*r 

xx  r 


♦  4r*»(t/i6  -  4/9t)  ♦  *r2[d/2  -  r(l  -  4/3*)]  2.  OP-1) 


Host  Investigators,  however,  will  neglect  the  loss  of  inertia  when  calculating  the 
bending  ?tr§ss  due  to  corner  radii  and  fsfyae  that  L  ■  bd  /.12,  and  the  error  in 
stress  becomes 


xx' r 


I 


-  If  100 


_  [  [btd-2r)s  +  2(b-2r)r3  ♦  6<b-2r)  (d-r)r2  ♦  48r**(^-  -  —•) 
*  M3  “ 


I2wr*  |  y  -  r 


-1.0  X  100 


(F-2) 


45  Chamfer 


Now  consider  Figure  Sb,  which  shows  a  rectangular  bean  with  45°  corner  chaafer 


c.  The  true  aouent  of  inertia  (I  )  about  the  x-x  axis  is: 

XX  c 

(IXX)C  -  (bd3/12)  -  (c2/9)[c2  ♦  (1/2) (3d- 2c)2]  , 
and  thus 


(F-3) 


-4c2 fc2  ♦  i  (3d- 2c) 2 

e  a  .  .  i  ■  -*T— .  * . • 


3bd3 


X  100. 


(F-4) 


The  errors  were  calculated  for  various  values  of  d/b  as  a  function  of  r/'d  fro® 
(F-2)  and  c/d  fro®  (F-4).  These  results  are  shown  in  Table  13. 


APPENDIX  6.  COMPUTER  ANALYSIS  WORKSHEET 


The  tables  in  this  report  should  suffice  to  permit  error  determinations.  A 
computer  program  is  available  at  MTL  to  expedite  such  computations  however.  NFL 
will  compute  the  error  analysis  upon  request  if  the  following  for®  is  filled  out 
as  completely  as  possible  and  nailed  to  Nr.  George  Quinn. 


Replies  will  be  kept  confidential. 


M-'  >  ,  >  i 


i. 

T 


—f  t— 


Plaasa  use  consistent  units  of  measure*  ive»»  all  dimensioni  in  inches*  millimeters* 
or  centiaeters,  ate. 

1.  Three-  or  four-point  flexure:  (circle  one) 

2.  Specimen  height  (thickness) 

3.  Spec  in en  width 

4.  Specimen  total  length  . 

5.  Specimen  edge  chamfer  radius  or  length*  (circle  one) 

6.  Fixture  outer  span 

7.  Fixture  inner  span  (if  4-point) 

8.  Fixture  outer  bearing  (s),  radius 

9.  Fixture  inner  bearing (s) ,  radius 


R 


10.  Precision  of  the  micrometer  used  to  measure 

specimen  height  and  width 

11.  Accuracy  of  the  fixture  spans 

12.  Specimen  twist  or  lack  of  parallelism  of  two  faces  $ 

13.  Fixture  twist  $  or  lack  of  parallelism  of  bearings 

14.  Length  of  bearing  fixture  O 

15.  Accuracy  of  centering  the  inner  bearing  (s)  relative 

to  the  outer  bearings 

16.  Do  your  fixture  bearings  rotate?  Or  are  they 

fixed*  such  as  knife  edges?  (circle  one) 

17.  Do  your  fixture  bearings  articulate  to  accommodate 

specimen  twist  or  warpage?  (circle  one) 


£ 2S8JJL 


yes 


18.  Error  in  measured  break  load 


H  Length 


Chamfer 


Specimen  Non- 
parallelism 
(Twist) 
Error 


fixed 


JSL 


Fixture  '<£>* 

Non- parallel Isa 
(TWIst)  Error 
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DEPARTMENT  OP  THE  ARMY 

UJ.  ARMY  LABORATORY  COMMANO 
MATIRIALS  TICHNOLOGY  LABORATORY 
WATIRTOWN,  MA88ACHUSITT8  02172-0001 


SLCMT-IMA-T  6  May  1988 

MEMORANDUM  FOR:  SEE  DISTRIBUTION 

SUBJECT:  Technical  Report  MTL  TR  87-35,  "Errors  Associated  with  Flexure 
Testing  of  Brittle  Materials" 

Errata  sheet  for  subject  report: 

Page  6,  Table  4: 

Change  E/a^  ■  500  to  E/a^  *  1000. 

Page  13,  Equation  18a: 

Add  parentheses  around  Ofc/2  (a^/2) . 

Page  14,  Equation  18b: 

Change  the  quantity  (e/o^)  to  (E/aj,) . 

Page  23,  Table  14: 

'  The  second  line  from  the  bottom  of  the  table;  change  p  to  P. 

Page  37,  Equation  C-14: 

Change  the  quantity  (b/J,')2  to  (n/i*)2. 

Page  41,  Equation  F-l: 

Change  the  quantity  ir/16  to  ir/18. 


FOR  THE  COMMANDER: 


•1  End .  ■> 
Diat-Liot 


DIANE  VALERI 
Chief 

Publishing  &  Visual  Information  Bi . 


